NPS ARCHIVE 
1998.09 
POLLARD, S. 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 




THESIS 



DEVELOPMENT AND VERIFICATION OF AN 
AERODYNAMIC MODEL FOR THE NPS FROG UAV 
USING THE CMARC PANEL CODE SOFTWARE SUITE 

by 

Stephen J. Pollard 
September 1998 

Thesis Advisor: Max F. Platzer 

Co-Advisor/Second Reader: Kevin D. Jones 



Thesis 

P6875 



Approved for public release; distribution is unlimited. 



DUDLEY KNOX LIBRARY 
MAVAL ROSTGRADUA1 F SC HO 
\’0-NTr=RE> CA ' A 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-01 o „ 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruction 
searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Senn 
comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, V 
Washington headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, V 
22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

September 1998 Aeronautical Engineer Thesis 


4. TITLE AND SUBTITLE 

DEVELOPMENT AND VERIFICATION OF AN AERODYNAMIC MODEL FOR THE NPS 
FROG UAV USING THE CMARC PANEL CODE SOFTWARE SUITE 


5. FUNDING NUMBERS 


6. AUTHOR(S) 
Pollard, Stephen J. 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 
Monterey, CA 93943-5000 


8. PERFORMING 
ORGANIZATION REPORT 
NUMBER 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 


10. SPONSORING/ 



MONITORING 

AGENCY REPORT NUMBt'K | 

11. SUPPLEMENTARY NOTES 

The views expressed in this thesis are those of the author and do not reflect the official policy or position of the Department of 
Defense or the U.S. Government. 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Approved for public release; distribution is unlimited 

13. ABSTRACT (maximum 200 words) 

The CMARC panel-code is evaluated for the development of an aerodynamic model of the Naval Postgradua e 
School FROG Unmanned Air Vehicle (UAV). CMARC is a personal computer hosted panel-code software suite for 
solving inviscid, incompressible flow over complex three-dimensional bodies. A panel model of the NPS FROG UAV 
is developed to obtain stability derivative data at the cruise flight condition. Emphasis is placed on comparing the 
CMARC data to aerodynamic models obtained from classical design techniques and parameter estimation. Linearized 
longitudinal and lateral-directional state-equation models are used to compare the dynamic response of each data set. 
In addition, CMARC is used to generate static-source and angle-of-attack sensor position corrections. Position ! 
corrections are provided in look-up table and curve-fit formats. The aerodynamic model obtained with CMARC i| 
demonstrated higher fidelity dynamic longitudinal response than the classical design model. Dynamic lateral - 
directional response is similar to that obtained from classical design techniques. Adjustment through comparison with 1 
flight-test data is still required to optimize the CMARC model. Future studies should concentrate on improving i 
CMARC modeling of fuselage side force through the addition of wake separation lines. Additionally, the propellu 
disk should be modeled in an attempt to capture the effects of increased dynamic pressure over the horizontal and 



vertical tail surfaces. 



14. SUBJECT TERMS 

Unmanned Aerial Vehicles, UAV, CMARC, Panel Method Ames Research Center, PMARC, Panel Code, 
Stability Derivatives, Boundary Layer Code, Aircraft Dynamic Response 


15. NUMBER O? 
PAGES 

195 


16. PRICE CODR 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION OF 
THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION 
OF ABSTRACT 

UL 



NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 



Prescribed by ANSI f-M 



12b. DISTRIBUTION CODE 



1 




11 



Approved for public release; distribution is unlimited 



DEVELOPMENT AND VERIFICATION OF AN AERODYNAMIC MODEL FOR 
THE NPS FROG UAV USING THE CMARC PANEL CODE SOFTWARE SUITE 

Stephen J. Pollard 
Commander, United States Navy 
M.S., Naval Postgraduate School, 1997 

Submitted in partial fulfillment of the 
requirements for the degree of 



AERONAUTICAL AND ASTRONAUTICAL ENGINEER 



from the 



NAVAL POSTGRADUATE SCHOOL 
September 1998 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHC' 

MONTEREY CA 330^3-5101 

ABSTRACT 



The CMARC panel-code is evaluated for the development of an aerodynamic 
model of the Naval Postgraduate School FROG Unmanned Air Vehicle (UAV). 
CMARC is a personal computer hosted panel-code software suite for solving inviscid, 
incompressible flow over complex three-dimensional bodies. A panel model of the NPS 
FROG UAV is developed to obtain stability derivative data at the cruise flight condition. 
Emphasis is placed on comparing the CMARC data to aerodynamic models obtained 
from classical design techniques and parameter estimation. Linearized longitudinal and 
lateral-directional state-equation models are used to compare the dynamic response of 
each data set. In addition, CMARC is used to generate static-source and angle-of-attack 
sensor position corrections. Position corrections are provided in look-up table and curve- 
fit formats. The aerodynamic model obtained with CMARC demonstrated higher fidelity 
dynamic longitudinal response than the classical design model. Dynamic lateral- 
directional response is similar to that obtained from classical design techniques. 
Adjustment through comparison with flight-test data is still required to optimize the 
CMARC model. Future studies should concentrate on improving CMARC modeling of 
fuselage side force through the addition of wake separation lines. Additionally, the 
propeller disk should be modeled in an attempt to capture the effects of increased 
dynamic pressure over the horizontal and vertical tail surfaces. 
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I. INTRODUCTION 



A. BACKGROUND 

Computational fluid dynamics (CFD) is increasingly used as a design and analysis 
tool. As the price of computer hardware drops and computational power increases, CFD 
becomes more attractive to a larger audience. CFD tools range from the high-end three- 
dimensional (3D) Navier-Stokes solvers for compressible, viscous fluids to potential flow 
solvers for incompressible, inviscid flows. This paper discusses the development of an 
aerodynamic model of the Naval Postgraduate School (NPS) Fiber Optic Guided (FROG) 
Unmanned Air Vehicle (UAV) using a personal computer hosted panel code. 

Validation of the Personal Simulation Works panel code software and initial 
FROG UAV modeling was completed in June 1997. Results are reported in Ref. [1], 
This study expands on the modeling effort to include the development of a complete 
aerodynamic model at the trim flight condition. Stability derivative data are obtained 
from the panel code and then compared to data from classical design calculations and 
parameter estimation. A linearized state space simulation is used to extract the air vehicle 
dynamic modes and to model control response. Results are compared for all three data 
sets. 

The Personal Simulation Works software suite, consisting of LOFTSMAN, 
CMARC and POSTMARC, is used for all aspects of the study. The software provides 
for panel model development, input file processing and the visualization of results. 
Emphasis is placed on verifying both the accuracy and suitability of the CFD programs 
for aerodynamic modeling. 

Until recently, personal computers (PC) did not have the computational power or 
memory to be practical for panel code CFD programs. Things have changed with the 
introduction of the Pentium class PC and low cost RAM. AeroLogic capitalized on the 
power of the Pentium class PC and developed Personal Simulation Works (PSW). PSW 
is based on the 3D, low order, inviscid potential flow solver named CMARC. CMARC is 
a re-hosted version of NASA’s Panel Method Ames Research Center (PMARC). 
PMARC was re-written in the C language and compiled for IBM compatible PCs. 
CMARC runs under the DOS operating system. It will also run in a DOS window under 
the WINDOWS 3.x, 95 or NT operating systems. Enhanced capabilities include; 
improved memory management, an expanded set of command line switches and 
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provisions for expanded boundary layer post-processing capabilities. However, the core 
processing algorithms remain the same as implemented in PMARC. 

LOFTSMAN, the PSW pre-processing program, is used to mesh complex 3D 
bodies and create input file patches. The program runs under the Windows operating 
system and allows the user to loft conics based 3D surfaces. The program automatically 
creates CMARC, PMARC or VS AERO input patches based on desired panel densities 
and distribution. 

POSTMARC is used for flow visualization and integration of resultant forces. It 
runs under the Windows 3.x, 95 and NT operating systems. POSTMARC reads CMARC 
or PMARC output files and provides for the visualization of model geometry, wake 
stepping, on and off-body streamlines and surface phenomena. 

B. REQUIREMENTS 

The Naval Postgraduate School Aeronautics Department is integrating UAV 
hardware and software to demonstrate autonomous flight, trajectory tracking and 
automatic landing. Closed-loop flight control development requires a valid aerodynamic 
truth model for the UAV airframe. The introduction of each new airframe requires the 
development of a new aerodynamic truth model. Most recently, Papageorgiou [Ref. 2] 
developed and tested an aerodynamic model for the NPS FROG UAV based on classical 
methods. His model produced a close match to flight test results in the longitudinal axis. 
However, the lateral-directional axis required modifications based on measured flight test 
data to produce acceptable results. With the availability of low cost panel code CFD 
capabilities, it is suggested that a panel code model of the FROG UAV will provide more 
refined stability derivative data for an initial aerodynamic model. 

In addition to a valid truth model, accurate pitot-static and angle-of-attack sensors 
are required for highly augmented flight control systems. CMARC is well suited for 
solving on-body static pressure distributions and off-body flow velocities over the 
predominately attached flow fields of a fuselage forebody. This proves particularly 
useful for generating pitot-static and angle of attack correction curves and look-up tables. 

C. STATEMENT OF OBJECTIVES 

The Naval Postgraduate School Aeronautics Department has both active CFD 
research and avionics development programs. The primary purpose of this investigation 



2 



is to verify the accuracy and suitability of the PSW software suite for developing an 
aerodynamic model of the NPS FROG UAV. Specific objectives are as follows: 

• Demonstrate panel code modeling, processing and visualization on a Pentium 
PC using the PSW software package. 

• Develop and analyze a panel code model for the NPS FROG UAV using PSW 
to obtain a complete set of stability derivatives at the cruise trim condition. 

• Verify the aerodynamic model through comparison to data obtained from 
classical design calculations and parameter estimation. 

• Demonstrate techniques for producing angle-of-attack vane and pitot-static 
correction curves. 

• Develop a user guide that outlines PSW panel code modeling for the 
extraction of stability derivatives. 



3 






























































4 



II. OVERVIEW OF PERSONAL SIMULATION WORKS 



A. GENERAL 

Personal Simulation Works is a PC based software suite that provides for the three 
primary CFD requirements; 3D modeling of an aircraft (LOFTSMAN), panel code flow 
solver (CMARC), and post-processing of the computed flow Field (POSTMARC). The 
software package contains three applications hosted on the IBM compatible personal 
computer. Each software program is discussed separately. 

B. LOFTSMAN 

LOFTSMAN is a Windows hosted 3D modeling tool that generates surface panel 
distributions for CMARC or PMARC input files. The program is based on conics, which 
allows rapid lofting of streamlined bodies such as aircraft fuselages and engine nacelles. 
In addition, wing and control surfaces can be designed with the extensive library of airfoil 
templates or with user specified coordinates. The software is well documented, including 
a tutorial, in the Personal Simulation Works User Guide [Ref. 3]. LOFTSMAN is 
primarily designed for creating new objects, but an existing airframe can be matched 
quite closely with just a detailed three-view drawing that includes frame cross sections. 
Appendix A outlines the development of the FROG UAV model. 

1. Streamlined Bodies 

LOFTSMAN functionality is divided into Body Objects and Wing Objects. In 
general, they remain separate unless the intersection between a wing and body is required. 

Body Objects are created using a family of curves called second-degree conics. 
Circles, ellipses, parabolas and hyperbolas are among this group. An entire fuselage is 
described by specifying just four lines. These are the top waterline (TW), bottom 
waterline (BW), the maximum breadth line (MB) and the waterline of the maximum 
breadth line (WW). For each line, the beginning, ending and a few points along the line 
are specified. Control points are also specified with a curvature factor that allows 
LOFTSMAN to generate a smooth conic between the points. The power of conic lofting 
will become evident when discussing the modeling of the complex FROG UAV fuselage 
in Chapter V. 
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Wings and Control Surfaces 



Wings and control surfaces are easily specified in LOFTSMAN using a short 
input file created with any text editor. The file specifies root, intermediate and tip rib 
section, location, axis, chord and incidence. LOFTSMAN then fairs a smooth surface 
through the rib sections. Washout is specified by varying the incidence of the root and tip 
ribs. Sweep-back is controlled by staggering the tip rib location with respect to the root 
rib. Once the general wing surface is specified, control surfaces such as ailerons, flaps 
and elevators can be deflected and meshed. 

3. Patches 

LOFTSMAN automatically meshes 3D surfaces and creates patches for 
CMARC/PMARC input files. The distinction between a mesh and a patch is important. 
A mesh is a set of quadrilateral and triangular panels that represent the surface of a wing 
or body. When the set of panels is organized and formatted to create a sub-component 
portion of a CMARC or PMARC input file, it is called a patch. 

A body or wing surface is first meshed at a density specified by the user. Panel 
compression options include cosine and half-cosine spacing. After meshing the object, 
one saves it to a text file as a formatted patch. One then opens the patch file with any text 
editor and copies/pastes the patch text into the appropriate location in the CMARC input 
file. 

Each control surface deflection requires a separate mesh and formatted patch. For 
instance, to evaluate roll performance one needs to separately mesh an upward aileron 
deflection on the right wing and a downward deflection on the left wing. If multiple 
deflections of a single control surface are required, each deflection must be meshed 
separately. 

C. CMARC 

CMARC is the C version of PMARC low-order, 3D panel code. Inviscid, 
irrotational, incompressible, potential flow is assumed. Low-order means that source and 
doublet strength distribution is constant across a panel. There is no attempt to match the 
source or doublet strength of an adjacent panel at a common edge. Advanced features 
include internal flow modeling and time stepping wake models. 
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PMARC version 12.19 was released as FORTRAN 77 source code in 1992. 
CMARC was rewritten in the C language and compiled for hosting on IBM compatible 
personal computers by AeroLogic, Inc. The program runs under the DOS operating 
system. It will also run in a DOS window from Windows 3.1, 95 or NT. Enhanced 
features include command-line options and flexible memory management. Command 
line options simplify batch processing by adding an extensive set of switches that can be 
set external to the CMARC input file. Flexible memory management provides for the 
automatic sizing of arrays without having to recompile the source code. 

D. POSTMARC 

POSTMARC is a Windows post-processing program for the visualization of 
CMARC and PMARC output files. Capabilities include body geometry, wake stepping, 
surface pressure and streamline visualization. POSTMARC also provides the capability 
to integrate pressure and skin friction forces over the model geometry. This proves 
particularly useful when one desires to recalculate loads around a different center of 
gravity. 

An interesting feature for design work is the integration of panel surface area to 
obtain total wetted area. After lofting a new geometry in LOFTSMAN, a quick check of 
geometry is made by running CMARC with the -g command line toggle. The total 
wetted area is then checked in POSTMARC. This function is particularly useful when 
working to reduce skin friction drag. 

Versions 1.17.3 and later of POSTMARC include the capability to integrate skin 
friction drag coefficient over the model geometry. It is important to note that a key piece 
of the drag equation is missing from a POSTMARC solution. CMARC provides induced 
drag from the surface pressure distribution and skin friction drag from the 2D boundary 
layer code. Skin friction is only calculated up to the point of boundary layer separation. 
Pressure drag due to separation, a major portion of the drag equation, is missing from a 
CMARC/POSTMARC solution. 

In fact, if one isn’t careful, POSTMARC drag calculations can be misleading. 
Take for instance two similar model configurations with only minor geometry differences 
that do not affect wetted area. It is possible for the model with more flow separation to 
have less skin friction drag because there is no CMARC output for skin friction 
coefficient after the boundary layer code predicts separation. During iterative design 



7 



work, this could lead to the incorrect conclusion that the design team is reducing overall 
drag. Perhaps a better function for LOFTSMAN than integrated skin friction drag would 
be a function that predicts the percentage of attached flow and laminar flow. Iterative 
design changes could be made that maximize laminar flow and minimize separated flow. 
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III. CMARC PANEL CODE THEORY 



A. POTENTIAL FLOW PANEL CODE THEORY (CMARC/PMARC) 

Potential flow theory involves the superposition of sources and doublets to 
generate the desired flow field around a 3D body. It assumes inviscid, irrotational and 
incompressible flow. As such, valid solutions are only obtained at low Mach numbers 
and for flow fields without large areas of separation. 

The basic concept of panel methods, as outlined by Bertin and Smith [Ref. 4], 
requires the modeling of the desired 3D configuration with a large number of 
quadrilateral and triangular panels representing the surface of the aircraft. A series of 
sources, doublets and vortices is then distributed on each panel. Superposition allows the 
simultaneous computation of the singularity strengths required to satisfy flow tangency 
on the surface. The inviscid, irrotational and incompressible flow field represented by the 
superposition of sources and doublets satisfies the Laplace equation: 



V 2 O = 0 



3.1 



Using Green’s Theorem, the potential at any point P in the flow is represented by: 






3.2 



Where (o-oj represents the potential from the doublet distribution and 
n ■ (v<t> - VO,.) represents the potential from the source distributions. 

CMARC is a low order panel code that assumes constant source and doublet 
strength distributions across each panel. Figure 3.1 shows a panel layout for a generic 3D 
wing fuselage configuration. It is important to note that for a 3D solution, there is an 
equivalence to surface doublet and surface vortex distributions. CMARC implements 
source and doublet distributions. 
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Figure 3.1 Typical Wing/Body Panel Code Configuration, from Ref. [5]. 

As mentioned previously, the general boundary condition imposed is tangential 
flow at the surface. CMARC, as outlined in Ref. [3], allows the modification of this 
boundary condition on individual panels or groups of panels. A normal surface velocity 
distribution may be specified to simulate flow into or out of ducts. 

In order to produce lift, a potential flow panel code requires a method to 
implement the Kutta condition. As noted in Anderson [Ref. 6], the Kutta condition at the 
trailing edge implies that the circulation, T, around an airfoil is such that the flow exits 
the trailing edge smoothly. In addition, the velocities leaving the top and bottom surfaces 
are finite and equal in magnitude and direction. 

Panel codes impose the Kutta condition by the shedding of wake panels along the 
trailing edges or separation lines. Wake panels are similar to a surface panel with only a 
doublet distribution. The doublet strength of the attached wake panel equals the 
difference in doublet strengths of the two adjacent surface panels. 
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The CMARC core panel code processing engine is functionally equivalent to the 
PMARC panel code module. The implemented equations are well documented by Ashby 
et al. [Ref. 5], The PMARC documentation includes a wing-body combination, shown in 
Figure 3.1, evaluated by PMARC with good correlation to experimental data. The 
results are shown in Figures 3.2 and 3.3. In addition, Lambert [Ref. 7] compared 
PMARC panel code results to several theoretical and experimental test cases with good 
correlation at low angle-of-attack. Sensitivity to wake placement is highlighted by his 
studies. 

Wake positioning can have a large influence on potential flow solutions. A wake 
is obviously attached to the trailing edge of wings and control surfaces with sharp, thin 
trailing edges to produce the Kutta condition. However, wake positioning on streamlined 
fuselages, missile airframes and nacelles is more of an art than science. Recently, Tuncer 
and Platzer [Ref 8] investigated generalized wake placement techniques for cylindrical 
bodies of revolution with good correlation to experimental data at up to 20 degrees angle- 
of-attack. The techniques were used with success by this author in Ref. [1] for the 
verification of CMAR.C calculations for flow over an inclined 6:1 prolate spheroid. 




Figure 3.2 Comparison of Experimental Data and PMARC Results for Two 
Spanwise Stations of the Wing/Body (a = 4°), from Ref. [5]. 
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Figure 3.3 Comparison of Experimental Data and PMARC Along the Fuselage 
Centerline of the Wing/Body Configuration (a = 4°), from Ref. [5]. 



B. CMARC BOUNDARY LAYER ANALYSIS THEORY 

CMARC and PMARC use the same two-dimensional integral method to calculate 
boundary layer characteristics along a surface streamline. A transition model 
automatically switches from laminar to turbulent calculations. The developers of the 
PMARC code chose a 2D integral routine over a 3D finite difference grid method due to 
speed and robustness of the calculations [Ref. 5]. Building a finite difference grid is a 
difficult and time consuming process requiring the user to develop grids over complex 3D 
surfaces. In addition, boundary layer calculation times can easily exceed that required for 
the basic potential flow solution. Reference [9] gives a good outline of three-dimensional 
finite difference methods. 



12 



The CMARC and PMARC User’s Guides [Refs. 3 and 5] contain detailed 
discussions on the development of the CMARC/PMARC boundary layer code starting 
from the two-dimensional momentum equation: 



dd , , 0 dU 1 

— + (2 +H)* = — C, 

dt] ’ U dr] 2 f 



3.3 

The momentum integral equation is numerically integrated along a surface streamline. 

The laminar region of the boundary layer is modeled by numerically integrating 
the following exact differential equation. The equation is solved iteratively through 
numerical integration along a streamline starting at a stagnation point [Ref. 5]: 




0.4 5 v 

U(i)‘ 



J(l + 2.222 g(K, u)p(n)" in + OiOY 

0 



r u(o) V 

yU(m\ 



3.4 



where: U - velocity at outer edge of boundary layer 

9 - momentum thickness 

K-^ 

v dr/ 

p - generalized coordinate along a streamline 

The value g(K,u) is based on exact solutions for a number of pressure 
distributions. Initial work was conducted by Thwaites with improvements by Curie 
[Ref. 10]: 

g(K,M) = F 0 {K) - juG 0 {K) - 0.45 + 6K 3.5 



CMAR.C uses an empirical transition model based on the average pressure 
gradient, K , for predicting laminar to turbulent transition. The following relations are 
used to calculate the transition point [Ref. 5]: 
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3.6 



\Kdrj 

K=— 

7" 7 

where rj ms is the streamline coordinate at instability. And, K is the local pressure 
gradient at boundary layer instability [Ref. 5]: 

K = -0.4709 + 0.1 1066* ln(Re*) + 0.0058591 * In 2 (Re* ) (o < Re* < 650) 

3.7 

K = 0.69412 - 0.23992 * ln(Re*) + 0.0205 * In 2 (Re*) (650 < Re* < 10000) 

The local Reynolds number at transition is correlated to K with the following expressions 
[Ref. 5]: 



K = -0.0925 + 0.00007* Re* 

K = -0.12571 + 0.0001 14286 * Re* 

K = 1.5938 1 - 0.45543 * ln(Re* ) + 0.032534 * In 2 (Re*) 



(0 < Re* < 750) 

(750 < Re* < 1 IOO) 3.8 
(l 1 00 < Re* < 3000) 



At transition, the initial turbulent shape factor, H, is given by the following empirical 
formula that is a fit to data developed by Coles [Ref. 10]: 



H = 



1.4754 

l°g,o( Re e) 



+ 0.9698 



3.9 



Provisions are made to check for turbulent reattachment if laminar separation is 
encountered. At laminar separation, a point calculation is made to determine if the 
boundary layer will reattach. If reattachment is predicted, the boundary layer code 
immediately switches to turbulent calculations. No attempt is made to model the laminar 
separation bubble or provide a transition length. After laminar separation is predicted, 
the following empirical relations are used to determine if reattachment occurs [Ref. 5]: 
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(Re, >125) 
(Re, <125) 



3.10 



K = 0.0227 - 0.0075 75 * Re, - 0.00000 1 1 5 7 * Re, 2 



K = -0.09 



The boundary layer code in CMARC uses a point transition model. No attempt is 
made to model a more representative transition length. Turbulent calculations begin at 
transition using the Nash-Hicks model [Ref. 5]. Calculations continue along the 
streamline until turbulent separation is predicted or the end of the streamline is reached. 
No boundary layer data is available after separation. 

The authors of PMARC caveat that their boundary layer calculations are quite 
accurate for predominately 2D flow but break down in regions of large cross flow near 
separation. This premise was tested in Reference [1] by comparing the predominately 2D 
flow over the inboard region of a high aspect ratio wing to the finite difference 
calculations performed by the Naval Postgraduate School Unsteady Potential Flow Code 
(UPOT). In general, CMARC provided correct trends for both the transition and 
separation points. However, in all cases, CMARC predicted early transition and late flow 
separation. As expected, the differences were greatest at lower Reynolds numbers where 
boundary layer thickness is larger. 

Reference [1] also modeled flow over an inclined prolate spheroid. The 6:1 
prolate spheroid was chosen because extensive experimental data is available. In 
addition, the three-dimensional flow around a prolate shperoid is similar to a streamlined 
slender fuselage. With proper wake placement, CMAJRC was found to produce accurate 
normal force and pitching moment coefficients. Over the three-dimensional body, 
CMARC boundary layer calculations also predicted early transition and late flow 
separation. Despite inaccuracies, CMARC boundary layer calculations remained useful 
when used as a design tool for visualizing the trend in transition and separation points 
with configuration changes. 
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IV. AERODYNAMIC MODEL OF THE FROG UAV 



A. BACKGROUND 

The Naval Postgraduate School Aeronautics Department is integrating UAV 
hardware and software to demonstrate autonomous flight, trajectory tracking and 
automatic landing. A core requirement for flight control law development is a valid 
aerodynamic truth model for the UAV airframe. A panel code model of the FROG UAV 
is one method for estimating many of the stability derivatives required for an 
aerodynamic truth model. This development effort concentrates on finding the static, rate 
damping and control power derivatives for both the longitudinal and lateral-directional 
axes. 

Panel code modeling usefulness goes beyond the development of aerodynamic 
coefficients. Flight control systems require accurate pitot-static and angle-of-attack 
sensor inputs. CMARC accurately solves on-body static pressure distributions and off- 
body flow velocities over the predominately attached flow fields of fuselage fore bodies. 
In this study, correction curves are generated for static source and angle-of-attack probe 
position errors. 



B. FROG UAV DESCRIPTION 

The FROG UAV is a small single engine flight test vehicle used for autonomous 
flight research by the Naval Postgraduate School Aeronautics Department. The aircraft 
was originally designated the FOG-R by the U. S. Army. It was designed as a small 
lightweight, battlefield observation platform that could be guided by a fiber optic data 
link. Table 4.1 presents the basic aircraft specifications. 

The aircraft is somewhat unconventional in that the engine is mounted in a nacelle 
tractor style above the fuselage and wing. The aft fuselage consists of a 1.75 inch 
diameter aluminum tube which connects the tail surfaces to the main fuselage. Figure 4. 1 
displays a three view drawing of the FROG UAV. 
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PARAMETER 


MEASUREMENT/UNITS 


Length 


8.125 ft 


97.5 in 


Height 


1.75 ft 


21 in 


Weight 


67.7 lbs 


Power Plant 


12 Hp/ 2 Cycle 


Wing Airfoil 


NACA 2415 


Horiz. Stab. Airfoil 


NACA 0006 (Approx.) 


S w (Sref) 


17.57 ft 2 


2530 in 2 


St 


3.174 ft 2 


457.1 in 2 


Sv 


0.9818 ft 2 


141.4 m 2 


C 


1.66 ft 


20 in 


C t 


0.958 ft 


11.5 in 


b w 


10.54 ft 


126.5 in 


b t 


3.313 ft 


39.75 in 


b v 


1.25 ft 


15.0 in 


It 


4.44 ft 


53.25 in 


lv 


4.44 ft 


53.25 in 


AR W 


6.32 


AR t 


3.46 


AR V 


1.59 


v H 


0.49 


V v 


0.16 



Table 4.1 FROG UAV Characteristics, after Ref. [2]. 
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Figure 4.1 FROG UAV Three-View Drawing. 
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The FROG UAV, as operated by NPS, is equipped with airspeed, angle-of-attack, 
altitude and control surface sensors. In addition, a miniature Inertial Measurement Unit 
(IMU) captures aircraft attitude, acceleration and body rates. Data is down linked to a 
mobile SGI workstation through a spread spectrum modem. Onboard GPS provides 
differential GPS navigation capability with the ground station used as a reference. The 
aircraft can be flown by conventional radio control or by up-linking flight control 
commands from the SGI workstation. 

Current flight control development requires a linearized aerodynamic model 
around the cruise trim point of 60 m.p.h. or 88 fit/s. This flight condition is selected for 
the development of stability derivative data with CMARC. Table 4.2 lists the aircraft 
parameters for the trim flight condition. 



PARAMETER 


MEASUREMENT 


UNITS 


Weight 


67.73 


lbs 


IXX 


12.52 


slug-ft 2 


IYY 


8.43 


slug-fit 2 


IZZ 


18.55 


slug-ft 2 


Airspeed 


60/88 


mph and ft/s 


Altitude 


800 


ft MSL 


Air Density 


0.002327 


slug/ft 3 


Center of Gravity 


34.5% 


M.A.C 


Cl trim 


0.4295 


n/a 


^ trim (est) 


-0.04 


degrees 


S Etnm 


5.1 


degrees 



Table 4.2 FROG UAV Trim Flight Condition, after Ref. [2]. 

C. FROG UAV MODELING 

1. General 

LOFTSMAN is utilized for the creation of all CMARC input file patches except 
for wing tips. In some cases, CMARC’s more efficient built-in capability to model 
standard NACA 4-digit wing surfaces could be used. However, with the requirement to 
mesh deflected control surfaces, all patches are created with LOFTSMAN from the start. 
The complete FROG UAV model with all patches activated is displayed in Figure 4.2. 



20 



Appendix A contains detailed descriptions of obtaining stability derivative data in a less 
formal user guide format. 




Figure 4.2 FROG UAV Panel Code Model with All Patches Activated. 

Some assumptions are made to simplify the modeling process. First, the 
horizontal and vertical stabilizers are modeled with a NACA 0006 section. The actual 
surfaces are constructed with a flat section, rounded at the nose and tapered starting at the 
control surface hinge line to a sharp trailing edge. The NACA 0006 provides a close 
approximation and allows the use of LOFTMAN’s built-in wing lofting capability. For a 
potential flow solution, this simplification is considered minor. 

A second simplification is made regarding the vertical stabilizer’s tip rib 
orientation. The actual rib is canted down 5° with respect to the longitudinal waterline. 
LOFTSMAN will only model a chord line that is parallel to the waterline (constant BL). 
Thus, for this study, the vertical tail tip rib is modeled with a constant BL, but the span is 
adjusted to maintain the same overall surface area. 

A third assumption is made regarding aileron control surface hinge lines. The 
actual FROG UAV ailerons are piano hinged along the upper surface. For simplification 
of the LOFTSMAN model, all surfaces are hinged at the centerline. 



21 



Finally, there is no attempt to model the tricycle landing gear struts or wheel 
assemblies. The landing gear components do not contribute significantly to the 
aerodynamic stability derivatives. However, they certainly need to be taken into account 
when measuring moments of inertia for a dynamic model. 

2. Modeling Coordinate System 

The model is developed using a coordinate system selected to simplify fuselage 
measurements. The +x-axis starts even with the nose and runs aft along the bottom of the 
fuselage, parallel with the tail boom. The bottom of the fuselage is used as the waterline 
with +z-axis in the up direction. This allows for easy vertical measurements when the 
aircraft is placed flat on a horizontal surface. The +y-axis runs from centerline outboard 
parallel to the right wing. Figure 4.3, which displays static-pressure source and alpha 
vane locations, also shows the location and origin of the modeling coordinate system. 

3. LOFTSMAN Patches 

LOFTSMAN is used to generate all the model patches except for wing tips. 
CMARC’s built-in capability is used to create wing tip patches. Appendix B contains 
listings of all the LOFTSMAN input files. Once a surface is meshed, the mesh is saved 
to a file as a CMARC/PMARC patch. The resulting text file is then opened, and the text 
is copied and pasted with any text editor into the patch definition section of the CMARC 
input file. LOFTSMAN patch files are not listed because they are redundant with the 
patches in the final CMARC input file listed in Appendix C. Multiple patches for the 
same control surface may be stored in the input file. For instance, two horizontal 
stabilizer patches with different elevator deflection angles may be kept in the CMARC 
input file. However, only one may reside in the active patch input section. All other 
patches for the same aircraft component must reside after the patch with the "last patch" 
toggle set (TNODS=5). 

When saving a patch, LOFTSMAN automatically takes care of all CMARC input 
formatting. A patch, as formatted by LOFTSMAN, assumes additional patches will 
follow in the CMARC input file. Therefore, the last segment’s TNODS variable is set 
TNODS=3. When the patch is the last patch in the input file, the TNODS variable must 
be manually set to TNODS=5. If CMARC hangs up while reading in geometry 
information, most likely TNODS=5 is missing on the last patch. 
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a. Fuselage Model 

The fuselage is lofted as a B-type body. A B-type body is used when 
major portions of the fuselage have a circular or oval cross section. The input file is 
listed in Appendix B. Only the right side is meshed, with a symmetric left side created 
by toggling the IPATSYM variable to IPATSYM=1. LOFTSMAN assumes that B-type 
bodies converge to a specific point at the fore and aft ends. The flat aft fuselage face does 
not provide this single point. A slight modification was made to the aft face to allow 
automatic meshing as a B-type body. The center of the aft face is extended very slightly, 
approximately 1/8 inch, to provide a convergence point for the final rear triangular 
panels. This small deviation is assumed not affect the aerodynamic fidelity of the model 
for a potential flow solution. 

The right side was originally meshed separately from the wing as a 20 x 
20 panel patch. This created a low order fit when the wing patch was butted to the side of 
the fuselage, resulting in overlapping panels. A final mesh was created that flowed 
around the wing root and fuselage intersection for a high order fit. All the fuselage 
panels at the wing root join with the adjacent wing panels. This mesh requires that the 
fuselage be broken up into six separate panels per side. They are the nose patch, the 
forward transition patch, the top and bottom wing root patches, the aft transition patch 
and finally the rear fuselage patch. Some manual editing is required to straighten out 
panels on the upper fuselage patch. When the six patches are added together, the final 
configuration is modeled with a 44x15 panel patch. 

b. Main Wing Patch 

The NACA 2415 wing is created with four separate patches to allow the 
addition of a deflected aileron mesh. CMARC comes with a broad selection of “*.SD” 
airfoil template files that are automatically loaded during installation. The 
“NACA2415.SD” file is used for this model. The 10x15 inboard wing patch runs from 
the wing root, past the flaps, to the start of the aileron. The 10x15 mid-wing patch covers 
the portion of the wing spanned by the aileron. The tapered outboard wing extension is 
made with a 6x15 patch. Finally, a 4x15 semi-circular wing tip patch is created using 
CMARC’s built-in wing tip functionality. The wing is set to a 4.5° incidence in the 
LOFTSMAN input file. Alternatively, the patch could be created with no incidence and 
then the patch coordinate system could be rotated in the CMARC input file. Together, 
the three wing patches add to make a 30x15 panel wing model. 
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The ailerons are meshed at zero degree deflection and five degrees of right 
rolling moment deflection. For the no deflection case, a single right wing patch is 
meshed with CMARC's symmetric patch toggle (IPATSYM=1) creating the left patch. 
For deflected ailerons, separate patches are meshed for each wing. The right wing aileron 
deflection is five degrees trailing edge up (T.E.U.) and the left wing five degrees trailing 
edge down (T.E.D.). 

c. Horizontal Stabilizer Patch 

The horizontal stabilizer and elevator patch is created with a single 10x16 
mesh using the “NACA0006.SD” airfoil template. As stated earlier, the actual UAV 
horizontal stabilizer has flat upper and lower surfaces and a rounded leading edge. The 
NACA 0006 airfoil is substituted because it provides a close approximation to the leading 
edged radius and thickness ratio. No other special modifications are required. A tip 
patch is not added because some of the resulting panels would be too small. In particular, 
the triangular panels closing out the aft end of the tip are too small in proportion to the 
other panels, creating singularities due to machine resolution. Initially, an attempt was 
made to model horizontal and vertical stabilizer wing tips, but the model would not 
converge with them. Leaving off tip patches does not significantly influence results 
according to the CMARC User's Guide [Ref. 3]. 

The elevator spans the entire horizontal stabilizer. Two patches are 
meshed, one with zero deflection and the second with positive five degrees (+5° T.E.D.) 
of deflection. The deflected patch is used to obtain the elevator control power derivatives 
(Clso CMSe and Cose)- 

d. Vertical Stabilizer Patch 

Two different versions of the vertical stabilizer are used for FROG 
modeling. An extended vertical stabilizer, which includes the "effective" area of the 
empennage tail boom, is used to measure the static sideslip and yaw rate damping 
derivatives. It is a 10x14 panel patch. The second patch models the actual area of the 
vertical stabilizer. It is a 10x12 patch and is meshed with a five-degree rudder deflection 
for measuring control power derivatives. 

The vertical stabilizer and rudder patches are created with a single mesh 
using the “NACA0006.SD” airfoil template. As with the horizontal tail stabilizer, the 
NACA 0006 airfoil closely approximates the vertical surface leading edge radius and 
thickness ratio. The LOFTSMAN input file is different in that a vertical wing surface 
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requires a modification to the rib axis. The rib axis must be specified with an x-axis 
rotation of 90°, a y-axis rotation of 0° and an unspecified (999.0) z-axis rotation. No 
symmetry is selected for the vertical stabilizer because the patch is already symmetric 
about the y=0 plane. As with the horizontal stabilizer, a tip patch is not added because 
some of the resulting panels would be too small. 

The rudder spans the entire vertical stabilizer. Two patches are meshed, 
one with zero deflection and the second with positive five degrees (+5° T.E.L.) of 
deflection. The deflected patch is used to obtain the elevator control power derivatives 
(Cls e, C\i5 e and CoSe)- 

e. Tail Boom Patch 

The tail boom patch is created as a single 12 x 10 mesh using a B-type 
body. Again, only the right side is meshed due to symmetry. The LOFTSMAN input file 
requires modifications at both ends in a similar fashion to the aft fuselage. A single point 
is added to allow convergence of the triangular panels at either end. With this point, the 
tail boom has the appearance of being tapered at both ends. The point is then manually 
edited out in the CMARC input file by replacing the “x” coordinate of the beginning and 
ending section panels with the correct value. In most cases, the tail boom is left out of 
solution to aid in convergence. This is due to the small overlapping panels at the fuselage 
tail boom junction. Being a slender, round tube directly in the fuselage slipstream, the 
tail boom is assumed to have a small influence on the stability derivatives. 

f. Engine Pod Patch 

The engine pod patch, or nacelle, is created as a single 15 x 10 mesh using 
a B-type body. Only the right side is meshed due to symmetry. The prop spinner is an 
integral part of the patch. No attempt is made to model the prop, engine heads or exhaust 
system. As with the tail boom, in most cases the engine pod is left off the model used to 
establish stability derivatives. 

g. Engine Pylon Patch 

The engine pylon patch is modeled with a single 15x10 mesh using an A- 
type body. A-type bodies are used to model surfaces similar to boat hulls with cornered 
surfaces or sharp chines. In addition, A-type bodies do not require the body to be 
completely enclosed. As a result, an A-body was selected to model just the sides of the 
pylon. Only the right side is meshed due to symmetry. A low order fit is achieved with 
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the adjacent fuselage and engine pod panels. This results in questionable pressure 
distributions. As a result, the pylon patch was turned off for most configurations. A 
future attempt could be made to create a high order fit between the other patches or to 
model the pylon as a flat surface. This will probably require manual editing of the 
intersecting patches. Even with a high order fit, a wake cannot be added to the trailing 
edge of the pylon. It would impact the vertical tail surface 

4. Common CMARC Input File Errors 

The patches created in LOFTSMAN are assembled into a single CMARC input 
file with any text editor. The default input file, which comes with CMARC, or any old 
input file may be modified. There are many errors that will cause CMARC to hang up 
without an error message. The two most common errors are forgetting to designate the 
last patch and incorrectly numbering the wake patches. 

The last patch must be designated by including a TNODS=5 setting in the last 
section of the last patch. If it is not included, CMARC hangs up when reading in the 
geometry. In a similar manner, the last wake must be designated with a NODEW=5 
setting. If the last wake is not designated, CMARC hangs up while reading in the wake 
information. Another common error involves incorrect wake to patch number 
association. Patch numbering changes whenever patches are disabled or reordered. The 
KWPACH field for each wake definition must be checked to make sure it reflects the 
current patch numbering. 

D. STATIC-PRESSURE SOURCE AND ALPHA VANE CORRECTIONS 

THROUGH OFF-BODY FLOW ANALYSIS 

CMARC is ideally suited for off-body flow analysis. Off-body streamlines may 
be placed through a point anywhere in the flow field. CMARC will then follow the 
streamline up and downstream the distance designated in the input file. This is 
particularly useful for flow visualization. In addition, CMARC calculates pressure 
coefficient and velocity at each point along the streamline. For this study, two 
streamlines are placed through the locations of the static -pressure source and alpha probe 
locations. Pressure coefficient is used to quantify static source position error and velocity 
is used to calculate alpha probe position error as a function of FROG UAV angle-of- 
attack. Both static pressure and AOA are digitized for down link to the ground station 
allowing the values to be easily corrected. Either a look-up table or curve fit correction 
can be applied subsequent to being passed to the flight control routines. This analysis 
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was previously published by the author in Ref. [1], but is included here to provide a 
single source of UAV modeling conducted to date. 

1. Description of the FROG UAV Pitot-Static and AOA Systems 

The pitot-static system and angle-of-attack probe share a common flight test 
boom extending from the nose of the UAV. The boom contains both the total and static 
pressure ports. Figure 4.3 depicts the general dimensions of the flight test boom 
installation and the modeling coordinate system. 

2. Modeling Off-Body Streamlines 

Streamlines are placed at the two locations indicated in Figure 4.3, which 
correspond to the static source and alpha probe locations. Two off body streamlines were 
activated in CMARC by setting NSTLIN=2 in the &SLIN1 line. Only a short distance of 
2 inches is selected up and downstream in the SU and SD fields to reduce the size of the 
output file. Figure 4.4 is a POSTMARC rendering of the two off-body streamlines used 
for sensor corrections. With the model at a t =0°, notice that the streamline is curving up 
at the angle-of-attack vane location 6.5 inches in front of the aircraft nose. 

3. Analysis of Static Source Position Errors 

The position error pressure coefficient, ACp pc or APp/q c , is a function of free 
stream Mach number and angle-of-attack provided that the static source is located outside 
of a thick boundary layer and sideslip is minimized [Ref. 11]. In the case of the FROG 
UAV with incompressible flow, APp/q c becomes a function of angle-of-attack only. 
Therefore, corrections can be simply defined as a function of measured angle-of-attack. 

A DOS batch file was executed to step the CMARC model through angles-of- 
attack ranging from -8° to 20°. The batch file incremented the angle-of-attack using 
CMARC ’s command line override feature. In addition, a new output file name was 
designated for each angle-of-attack. Position error pressure coefficient is then read from 
the off-body streamline listing of the output file at the location corresponding to the static 
source. Table 4.3 lists the values of APp/q c calculated from CMARC data. Figure 4.5 
displays ACp pc as a function of indicated angle-of-attack. The second order influence of 
angle-of-attack is clear with a second order curve fitting tightly through the data points. 
Of note, the error is relatively constant for a ±8° band around trim angle-of-attack. For 
incompressible flow, position error pressure coefficient is independent of airspeed and 
altitude. 
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Figure 4.4 FROG Off-Body Streamline Visualization with POSTMARC (a=10°). 

Position error pressure coefficient can be turned into position corrections for 
airspeed and altitude. The following relations were developed which assume small errors 
and incompressible flow: 

ViACp 

AV pc =— y 2 - and AV pc =V c -V i 4.1 

AV pc Vi 

AH pc =— and AH pc =H c -Hi 4.2 

^stdSo 

Where: 



AH pc is the altitude position correction. 

AV pc is the velocity position correction. 

AP d 

AC d = — — or position error pressure coefficient. 

q c 

a st( j is standard day density ratio. 
g 0 is the gravitational constant. 
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Table 4.3 displays corrections calculated for both airspeed and altitude at the 
FROG UAV trim condition of 88 ft/s and 800 ft MSL. The corrections are added to the 
indicated value to obtain the corrected value. Figures 4.6 and 4.7 display the corrections 
as a function of indicated angle-of- attack. Again, a second order curve fits nicely through 
the data points. Equations 4.1 and 4.2 can be used to implement a correction algorithm 
based on airspeed and altitude. 

4. Analysis of Alpha Vane Position Error 

Local flow field velocity is extracted from the off-body streamline listing to 
obtain local angle-of- attack. The alpha vane is assumed to capture the x-z component of 
the local velocity field and ignore cross flow in the y direction. Flow field velocity is 
turned into indicated angle-of-attack and angle-of-attack position correction with the 
following equations: 



a i 




180 

n 



degrees 



4.3 



o 

Aa pc =a t -ot{ degrees 



4.4 



A DOS batch file is executed to step the CMARC model, with an off-body 
streamline located at the vane position, through angles-of-attack ranging from -8° to 20°. 
Local velocity components are then read from the location corresponding to the alpha 
vane. Table 4.4 lists the values of Aa pc calculated from CMARC data. Figure 4.8 
displays Aa pc as a function of indicated angle-of-attack. Linear and second order curve 
fit equations are also indicated on Figure 4.8. Angle-of-attack correction is fairly linear 
through the FROG operating envelope, with approximately -1.25 degrees of position 
error at the FROG cruise trim condition. The corrections apply at all incompressible 
airspeeds and all altitudes. 

5. Summary of Off-Body Flow Field Analysis 

CMARC proved useful for both static-pressure source and alpha vane position 
corrections. Measured data may be corrected using look-up tables with the values in 
Table 4.3 and 4.4 or by using the curve fits in Figures 4.5 through 4.8. Flight testing is 
recommended for validation of sensor corrections obtained from this CMARC off-body 
flow field analysis. 
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Table 4.3 



UAV AOA 

a T (deg) 


ACp pc 

AP/q c 


V Correction 

AV pc =V c -Vj (ft/s) 


P Correction 

AH pc =H c -Hi (ft) 


-8 


0.1092 


4.8 


13.5 


-6 


0.1120 


4.9 


13.8 


-3 


0.1141 


5.0 


14.1 


-2 


0.1140 


5.0 


14.1 


-1 


0.1137 


5.0 


14.1 


0 


0.1132 


5.0 


14.0 


1 


0.1123 


5.0 


13.9 


2 


0.1111 


4.9 


13.7 


3 


0.1096 


4.8 


13.5 


4 


0.1078 


4.8 


13.3 


5 


0.1057 


4.7 


13.1 


6 


0.1034 


4.6 


12.8 


8 


0.0977 


4.3 


12.1 


10 


0.0909 


4.0 


11.2 


12 


0.0831 


3.7 


10.3 


14 


0.0741 


3.3 


9.2 


16 


0.0641 


2.8 


7.9 


18 


0.0530 


2.3 


6.6 


20 


0.0410 


1.8 


5.1 



Position Error Corrections for the NPS FROG UAV at V=88 ft/s and 
H=800 ft MSL. Derived from CMARC Panel Code Off-Body Flow 
Field Analysis. 
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ACp pc 



Source: CMARC Panel Code 
Trim Condition: V=88 ft/s H=800 ft MSL 




AOA Indicated (deg) 



Figure 4.5 Position Error Pressure Coefficient, AC P pw for the NPS FROG UAV. 

Derived from CMARC Panel Code Off-Body Flow Field Analysis. 
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AHpc (ft) 



Source: CMARC Panel Code 
Trim Condition V=88tt/s H=800 ft MSL 
Flight Test Boom Static Source (8.25" forward of nose) 




AOA Indicated (deg) 



Figure 4.6 Altitude Position Error, AH pc , for the NPS FROG UAV at V=88 ft/s 
and H=800 ft MSL. Derived from CMARC Panel Code Off-Body 
Flow Field Analysis. 
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Source: CMARC Panel Code 
Trim Condition: V=88 ft/s H=800 ft MSL 
Flight Test Boom Static Source (8.25" forward of nose) 




AOA Indicated (deg) 



Figure 4.7 Airspeed Position Error, AVpc, for the NPS FROG UAV at V=88 ft/s 
and H=800 ft MSL. Derived from CMARC Panel Code Off-Body 
Flow Field Analysis. 
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UAV AOA 

a T (deg) 


Velocity at Alpha Vane 


AOA Correction 
Aa=a T -a 1 (deg) 


AOA Indicated 
ctj (deg) 


V x (ft/s) 


Vy (ft/s) 


V z (ft/s) 


-8 


80.92 


1.66 


-12.23 


0.60 


-8.60 


-6 


81.27 


1.65 


-8.65 


0.08 


-6.08 


-3 


81.60 


1.64 


-3.21 


-0.75 


-2.25 


-2 


81.67 


1.63 


-1.47 


-0.97 


-1.03 


-1 


81.71 


1.63 


0.28 


-1.20 


0.20 


0 


81.73 


1.62 


2.13 


-1.49 


1.49 


1 


81.73 


1.61 


3.93 


-1.75 


2.75 


2 


81.70 


1.60 


5.72 


-2.00 


4.00 


3 


81.66 


1.59 


7.51 


-2.25 


5.25 


4 


81.58 


1.58 


9.30 


-2.50 


6.50 


5 


81.48 


1.57 


11.08 


-2.75 


7.75 


6 


81.37 


1.56 


12.88 


-2.99 


8.99 


8 


81.07 


1.53 


16.43 


-3.46 


11.46 


10 


80.67 


1.51 


19.98 


-3.91 


13.91 


12 


80.17 


1.48 


23.50 


-4.34 


16.34 


14 


79.61 


1.46 


26.99 


-4.73 


18.73 


16 


78.93 


1.43 


30.47 


-5.11 


21.11 


18 


78.18 


1.39 


33.90 


-5.44 


23.44 


20 


77.34 


1.36 


37.31 


-5.75 


25.75 



Table 4.4 Angle-of Attack Vane Position Error Corrections for the NPS FROG 
UAV. Derived from CMARC Panel Code Off-Body Flow Field 
Analysis. 
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Source: CMARC Panel Code 
Trim Condition: V=88 ft/s H=800 ft MSL 
Flight Test Boom Alpha Vane (6.5" forward of nose) 




Figure 4.8 Angle-of- Attack Vane Position Error, Aa pc , for the NPS FROG UAV. 

Derived from CMARC Panel Code Off-Body Flow Field Analysis. 
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E. 



DEVELOPMENT OF STABILITY DERIVATIVES 



In this section, CMARC is used to develop the longitudinal and lateral-directional 
stability derivatives for the FROG UAV. The development effort focuses on finding the 
static, rate-damping and control-power derivatives. The results obtained from CMARC 
are compared to data sets obtained from empirical estimation techniques and flight-test 
parameter estimation. In the next section, the stability derivatives are entered into a 
dynamic model to find modal frequency, damping and response to control deflections. 

The potential flow analysis performed by CMARC does not provide accurate 
viscous drag values. Therefore, total drag is estimated from the flight-test power-off 
glide ratio and cruise thrust required. The User Guide developed in Appendix A 
describes in detail how to gather stability derivative data with CMARC. An abbreviated 
discussion is presented below. 

CMARC contains built-in functionality to integrate forces and moments in all 
axes over the surface of a body. Force and moment coefficients are automatically 
nondimensionalized based on the mean aerodynamic chord, reference wing area, semi- 
span and center of gravity location in the CMARC BINP9 input line. Coefficients are 
presented in both wind and body axes. Of note, CMARC uses the semi-span to 
nondimensionalize rolling and yawing moment coefficients. Most texts on stability and 
control, including Roskarn [Ref. 12] and Etkin [Ref. 13], nondimensionalize rolling and 
yawing moments by span. Rolling and yawing moments are nondimensionalized by span 
in this study. Table 4.5 summarizes the factors for normalizing moments and angular 
rates. All rolling and yawing moment coefficients presented in this study have been 
normalized with span by dividing the CMARC output by a factor of two. Table 4.5 also 
indicates the characteristic time, t*, employed for angle rate data reduction. 

In addition to differences in nondimensionalizing moments, CMARC uses the 
typical CFD axis system shown in Figure 4.9. For this study, all work is performed in the 
stability axis system. Figure 4.9 also illustrates the standard stability axes implemented 
in this study. The sign of CMARC roll and yaw moments need to be reversed. The 
direction for positive control deflections is also shown in Figure 4.9. All control surfaces 
are patched with positive deflections using the convention in Figure 4.9. 
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A potential flow solution will not produce satisfactory results for bodies with 
significant areas of flow separation. Therefore, CMARC models must be analyzed in the 
linear slope regions for valid results. 
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Note: 1) CMARC nondimensionalizes roll and yaw coefficients with b/2. 



Table 4.5 Nondimensionalized Moment and Rate Equations. 

For the static stability derivatives, the CMARC model is run at two different 
angles-of- attack and one sideslip angle with zero control surface deflections. The slopes 
of the force and moment coefficients are then taken to produce the Cl« and C ma 
longitudinal derivatives and the Cyp, Qp and C n p lateral-directional derivatives. 

For the control-power derivatives, the model is run at the trim condition with 
successive control deflections. The difference between the results with and without 
control surface deflections yield the C[_s e , CM5e and Cose longitudinal and the Cysr, Qsr, 
C n 6r, Qsa and C n sa lateral-directional control-power derivatives. 

Development of the damping derivatives is not as straight forward. The static 
derivatives are obtained with motion disabled. For the longitudinal damping derivatives, 
the model is run with oscillating vertical plunging motion to obtain the Cl and Cm a-dot 
terms. The lift and pitching moment coefficients are broken into real (in phase with 
AOA) and imaginary (out of phase with AOA) components. The imaginary components 
are due to a-dot effects. Next, the model is run with oscillating pitch motion to obtain 
the combined a-dot and pitch rate terms. Subtracting the a-dot influence obtained from 
the plunging motion isolates the pitch rate-damping term from the pitch motion. 
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CMARC Axis System 




Stability Axis System 



Figure 4.9 CMARC CFD Axis System Compared To Stability Axis System. From 
References [3] and [14] Respectively. 
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For the lateral-directional analysis, the (3-dot terms are assumed to be small. This 
allows the model to be run with just oscillating roll and yaw motion. As with the 
longitudinal test case, the imaginary or out of phase component yields the combined 13- 
dot and rate-damping data. With the (3-dot terms assumed to be small, the oscillating 
motion yields the Cy, Ci, and C n roll and yaw rate terms directly. 

Throughout this chapter, stability derivatives obtained from CMARC are 
compared to three other data sets. The first comes from classical analysis through the 
methods presented by Roskam in References [12] and [15]. The second set comes from a 
classical aerodynamic model developed by Papageorgio [Ref. 2] and further refined 
through parameter estimation by Engdahl. The resulting parameter estimation model 
provides a good match to observed flight characteristics and is currently in use by the 
NPS Department for Aeronautical Engineering for closed-loop flight control 
development. The third set of data is for the Cessna 1 72 at the cruise flight condition. 
The data is only presented for an order of magnitude comparison. 

1. Longitudinal Stability Derivatives 



two runs of the CMARC model. The model is first analyzed at an angle-of- attack 
corresponding to the estimated trim condition. In this case, a t =0° is selected for the first 
run. A second CMARC run is conducted with angle-of-attack incremented by one or two 
degrees. Cl and C m are then extracted manually from the data files. The slope of Cl and 
C m versus angle-of-attack provide the Clcx and C ma longitudinal derivatives. For this 
study, several angles-of-attack were analyzed to check consistency of the slope. In 
addition, cttrim is calculated from the lift curve slope and trim lift coefficient. For the 
longitudinal analysis, only half the model is analyzed. The symmetric calculation mode 
is selected by setting both RSYM=0.0 and IPATSYM=0 in the CMARC input file. 
Equations 4.5 through 4.7 are used for these calculations: 



Static Longitudinal Stability Derivatives 



Three basic longitudinal stability derivatives can be measured with just 
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Two FROG UAV model configurations are analyzed in a build-up 
approach to check results against classical calculations and flight-test data. Figure 4.10 
shows the CM ARC models. First, just the wing and horizontal tail are considered. The 
patches for all other surfaces and wakes are turned off and the wing root is extended to 
centerline. Then, the blended wing/fuselage and horizontal tail are analyzed. Values of 
Ci_ a and C ma for these two configurations are presented in Table 4.6. 

Classical design calculations are also performed to estimate C ma for 
comparison to CMARC results. Equation 4.8 is used for the calculation of Cmc^ 



Cm a — a w 



(h-had-Vn-i- 1- 



w 



_ds_ 

da/ 



4.8 



In classical design, the horizontal tail downwash derivative, ds/da, is generally selected 
from empirical data. Using a taper ratio of TR=1:1 and an aspect ratio of AR=6, 
ds/da=0.4 is selected from empirical charts in Ref. [16] for the FROG UAV 
configuration. Also, calculations for de/da=0.25 are included to illustrate pitching 
moment sensitivity to the downwash derivative. Table 4.6 lists Cmafor each 
configuration. 

The final static longitudinal parameter required is total aircraft drag. Drag 
coefficient plays an important role in long-period aircraft dynamics. Unfortunately, 
potential flow panel codes such as CMARC do not provide accurate total drag estimates. 
They can provide good induced drag predictions. And, if equipped with a boundary layer 
code like that contained in CMARC, they can provide integrated skin friction results. 
However, a large total drag contribution in the form of pressure drag due to flow 
separation is not accounted for. Total drag estimates are made below using the two 
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simple techniques shown in Equations 4.9 to 4.12. The first method is based on flight- 
test glide ratio. The second is based on cruise power required and estimated prop 
efficiency. Note that the selected prop efficiency is relatively low due to the small 
propeller diameter, high RPM and pusher configuration. The two methods provide drag 
predictions within 10% of each other. The results are averaged to Cd=0.065 and included 
in Table 4.6. 

Method 1: Lift-to-Drag Ratio (L/D=7 from flight-test) 




4.9 



D 



9.67 lbs 



= 0.0611 



qS 0.5*0.002327 lb-s 2 /ft 4 *S 8 2 ft 2 / s 2 *17.57 ft 



4.10 



Method 2: Cruise Power Setting (HP=5, r|p=0.35) 




D T r HP r * 550 ftlbs/s/ HP*q p /V 
qS qS qS 



qS 



4.11 



5 HP * 550 f-lbs/s/HP*035m ft/s 



4.12 



0.5*0.002327 lb-s 2 /ft 4 * 88 2 ft 2 !s 2 *\1.51 ft 2 



Average: Coave =(0.061 1+0.069)/2 = 0.065 
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Simplified Wing and Horizontal Tail 





Blended Wing and Fuselage with Horizontal Tail 



Figure 4.10 CMARC Models of the FROG UAV for the Determination of Static 
Longitudinal Stability Derivatives. 
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Flight-test data for the short-period and phugoid modes were used for 
longitudinal parameter estimation. Values for Cl<x and C m a based on preliminary 
parameter estimation work by Engdahl are presented in Table 4.6. 



METHOD 


CONFIGURATION 1 


STATIC 

Antrim 

laeg) 


I LONGIT 
C|-a 

(per raa) 


UDINAL 

^ma 
(per raa) 


PARAME 

C D 


TERS 

Cd« 

(per raa; 


CMARC 


Wing/Horiz Tail 


-0.87 


4.86 


-0.835 


n/a 


n/a 


Panel Code 


Wing/Fuselage/Horiz Tail 


-0.01 


4.85 


-0.413 


n/a 


0.266 


Classical 


Wing/Horiz Tail - 8s/6a=0.25 


-0.81 


4.85 


-1.00 


n/a 


n/a 


Design 3 


Wing/Horiz Tail - 5e/8a=0.40 


-0.82 


4.82 


-0.70 


n/a 


0.253 


Parameter 

Estimation 4 


Flying Aircraft 


n/a 


4.09 


-0.42 


0.0655 


n/a 




C-172 6 


n/a 


0.31 


-0.89 


0.31 


0.13 



NOTES: 1) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Zero lift wing incidence is +6.5° from the longitudinal reference line. 

3) Classical design after Ref. [12]. 

4) Parameter estimation from NPS flight test data by Engdahl. 

5) Average drag obtained from L/D and cruise power analysis from flight test data. 

6) C-172 data from Ref. [12]. 

Table 4.6 Comparison of FROG UAV Static Longitudinal Stability Derivatives. 

Good correlation of lift curve slope (Clo) to classical design is 
demonstrated by the two CMARC configurations in Table 4.6. Both techniques over 
estimated Clq compared to flight-test data. Neither the classical design technique nor the 
CMARC model accounts for control surface leakage or separation effects. Modeling 
actual flight control gap geometry may provide a closer match to experimental data. 

Excellent correlation is achieved for C m a between the blended 
wing/fuselage CMARC model and flight-test results. The low value of C ma implies that 
ds/da is greater than the 0.4 value selected from empirical charts in Reference [16]. The 
simplified wing/tail model resulted in a C ma that closely matches classical design 
predictions with ds/da=0.35. Clearly, CMARC modeling provides the value-added 
capability to capture fuselage influences that are missed with classical design methods. 
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An investigation was performed to see if improved wake definition could 
better capture the ds/da downwash derivative. Initially a flexible wake was selected. 
However, at the cruise angle-of-attack, the flexible wake impacted the horizontal tail 
causing inaccurate results. A further investigation was made into a streamline based 
wake definition. A study by Walden et al. [Ref. 17] studied wake turbulence by 
modeling an aircraft flying in trail of a wake-generating wing. A horizontal tail trailing 
the main wing is a similar configuration. The study found that a streamline-based wake 
is the best method for modeling downwash effects. In the study, the model was initially 
run with a rigid-wake. A streamline was also defined downstream of the main wing 
trailing edge. After analyzing the first results, the wake was then predefined to follow the 
streamline. The result was a first iteration on modeling a flexible wake with a predefined 
rigid-wake. This technique was tried during the course of FROG UAV modeling. 
However, FROG geometry proved unsuited to this technique. The streamline trailing the 
wing ended up impacting the horizontal tail. This indicated that a predefined wake would 
also impact the horizontal tail, resulting in inaccurate result. Figure 4.11 shows the 
streamline impacting the horizontal tail at cruise angle-of-attack with rigid-wakes 
selected. With the already close correlation of the rigid-wake results to flight-test data, 
further efforts to investigate C ma were abandoned. 




z 

1 

y x 



Figure 4.11 Rigid-Wake Streamline for Investigation of Streamline Based Wake 
Definition. 
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In summary, for the static longitudinal stability derivatives, CMARC 
produced accurate values for C ma and a slightly high value of Ci_ a . Drag coefficient is 
obtained by averaging the values from the lift-to-drag ratio and cruise thrust required 
techniques. 



b. Longitudinal Damping Stability Derivatives 

The two aircraft motions illustrated in Figure 4.12 are used to develop the 
longitudinal rate-damping derivatives. A sinusoidal vertical plunging motion isolates the 
a-dot effects. Oscillatory motion is controlled with the CMARC BINP8A input file line. 
The oscillating pitch angle captures both the a-dot and pitch rate-damping terms. It is 
controlled with the CMARC BINP8B input file line. The a-dot influence developed 
from the plunging motion is subtracted from the pitching motion to isolate the pitch rate 
effect. All motion is conducted at a frequency of 2n rad/s, which equates to a reduced 
frequency of k=0.0595 for this configuration and trim airspeed. 




Figure 4.12 Aircraft Motion for the Determination of a-dot and Pitch Rate 
Damping. 
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The two models shown in Figure 4.13 are used for the rate-damping 
investigation. The complete wing/fuselage horizontal tail model is used for both the a- 
dot and pitch rate-damping terms, while the simplified horizontal tail only model is used 
to obtain a first approximation of pitch rate-damping. The User Guide in Appendix A 
provides a more detailed description of the techniques used to gather the dynamic 
longitudinal stability derivatives. 




Simplified Horizontal Tail Model 




z x 

Blended Wing and Fuselage with Horizontal Tail 

Figure 4.13 CMARC Models of the FROG UAV for the Determination of 
Longitudinal Rate-Damping. 
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Traditional design calculations, from Roskam [Ref. 12], are used for 
comparison to CMARC and flight-test results: 



r =2*C *n *V 

^ ^ La H IH r H 



ds 

da 



4.13 
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= -2*C Ln *t, H 



*V -L* 

H — 

c da 



4.14 



r =2*C *n * V 

L q L ^La H IH V H 

r = -J * r *n *V — 

where the dynamic pressure ratio at the tail is taken as r| H =l .0 

The dynamic longitudinal stability derivatives are presented in Table 4.7. 
CMARC-obtained values are compared to classical design values and C-172 stability 
derivatives. Flight-test data is unavailable for the damping derivatives. The aerodynamic 
model currently in use by the FROG research team implements longitudinal damping 
coefficients obtained from classical design techniques. The CLa-dot and C mq obtained 
from CMARC for the complete model provided a good match to the classical design data. 
However, both CLa-dot and C mq are over estimated by approximately 50% as compared 
to classical design values. One possible explanation for the additional damping from 
CMARC is that the classical approximation only estimates the horizontal tail influence on 
damping. CMARC also measures wing and fuselage damping influence. The horizontal 
tail by itself produces values close to the classical calculations. This supports the 
assertion that CMARC also captures the wing and fuselage influence for the whole 
model. 



4.15 



4.16 
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METHOD 


CONFIGURATION 1 


DYNAMIC LONGITUDINAL PARAMETERS 


C La-dot 

(aeg) 


C Lq 

(per raa; 


c 

^ma-dot 
(per raa) 


Cmq 
(per raaj 


Panel Code 


Blended Wing-Fuselage/Horiz Tail 


1.42 


6.82 


-6.24 


-11.78 


Horizontal Tail only 


n/a 


4.37 


n/a 


-11.94 


Classical 2 


Wing / Horizontal Tail 


1.56 


3.89 


-4.14 


-11.39 


Par. Est 3 


Flying Aircraft 


n/a 


n/a 


n/a 


n/a 




C-172 4 


1.7 


3.9 


-5.2 


-12.4 



NOTES: 1 ) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Classical design after Ref. [12]. 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref. [12]. 

Table 4.7 Comparison of FROG UAV Dynamic Longitudinal Stability Derivatives. 



c. Longitudinal Control-Power Derivatives 



The elevator control-power derivatives are obtained by substituting a 0° 
deflection horizontal-tail patch for one with +5° trailing edge down deflection. Only one 
run is required. The difference between the trim condition and the deflected value is 
divided by the elevator deflection as shown below: 
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Elevator control-power derivatives are presented in Table 4.8. A classical 
design estimate of elevator control-power is also shown for comparison. The following 
relationships from Roskam [Ref. 15] are used for estimating elevator control-power: 
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4.20 
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L Sc L 8p g 
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4.21 



where Clsf is the variation of lift coefficient with elevator deflection found from charts in 
Roskam [Ref. 15]. 

In general, the elevator control-power derivatives obtained with CMARC 
correlate well with classical design techniques. For both Cue and Cms e , CMARC derived 
elevator power is approximately 15 percent higher than obtained from classical 
calculations. However, CMARC significantly under-estimates elevator control-power 
compared to the parameter estimation model. The most likely explanation for this 
phenomenon is that the CMARC solution does not model the effects of prop wash. The 
higher dynamic pressure over the horizontal stabilizer from prop wash increases elevator 
effectiveness. The parameter estimation model captures this phenomenon. This same 
trend can be seen in the modeling of rudder control-power. An attempt could be made to 
model the prop disk effects using CMARC's normal panel velocity definition capability. 



METHOD 


CONFIGURATION 1 


ELEVATOR CONTRO 


POWER 


^Loe 

(per rad) 


Co6e 

(per rad) 


^mSe 

(per rad) 


Panel Code 


Blended Wing-Fuselage/Horiz Tail 


0.438 


0.01 


-1.199 


Classical 2 


Wing / Horizontal Tail 


0.39 


n/a 


-1.04 


Par. Est 3 


Flying Aircraft 


1.13 


n/a 


-1.62 




C-172 4 


0.43 


0.06 


-1.28 



NOTES: 1) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Classical design after Ref. [12]. 

3) Parameter estimation from NPS flight test data by Engdahl. 

4) C-172 data from Ref. [12].. 

Table 4.8 Comparison of FROG UAV Elevator Control-Power Derivatives. 
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2. Lateral Directional Stability Derivatives 



a. Static Lateral-Directional Stability Derivatives 

Development of the static lateral-directional stability derivatives is not as 
straightforward. Both sides of the airframe must be modeled by setting both RSYM-1.0 
and IPATSYM=1. This creates symmetric patches around the y=0 plane allowing 
CMARC to perform asymmetric calculations around the entire body and significantly 
increases processing times. 

Initially, a lateral-directional solution was attempted with a complete 
airframe model. The longitudinal model was modified by adding the vertical stabilizer. 
This required what was perceived to be a minor wake modification. Figure 4.14 
illustrates the configuration. The wing wake was terminated at the wing root to prevent 
interference with the vertical stabilizer. This left a wake gap for the vertical stabilizer. 
The engine nacelle and pylon patches were turned off because the wakes would also 
impact the vertical stabilizer. This configuration met with mixed success. The values of 
both Cyb and C n b were opposite the expected directions. Under closer inspection, it was 
determined that the false vortex shed from the wing root caused a destabilizing influence 
on the vertical stabilizer. Figure 4.15 illustrates the vortex using CMARC's built in off- 
body streamline capability. 

In order to obtain satisfactory results, the model is broken up into two 
separate groups. The aircraft is modeled with the blended wing and fuselage as one 
group and the horizontal and vertical stabilizers as another group. Separate solutions are 
obtained for each group and summed through superposition to obtain the whole aircraft 
solution. The two separate models are shown in Figure 4.16. This method should be 
used as a last resort. A complete solution should be used if the airframe geometry allows 
adequate wake separation from the vertical stabilizer. A complete solution would capture 
fuselage, wing and vertical stabilizer interactions. 

Once the model is defined, it is checked for lateral-directional balance at 
zero sideslip angle (yaw angle=0°). The side force, rolling and yawing coefficients 
should be zero when a trial run is performed at zero sideslip. If lateral-directional forces 
or moments are present, the model and wake geometry should be checked for symmetry. 
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Figure 4.14 Unsuccessful Lateral-Directional Model of the FROG UAV. False 
Wing Root Vortex Caused Destabilizing Influence on the Vertical Stabilizer. 




Figure 4.15 Destabilizing Wing Root Vortex Visualized with Off-body Streamlines. 
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Figure 4.16 Final Simplified Lateral-Directional Models of the FROG UAV. 
Separate Solutions are Summed to Obtain Complete Airframe Solution. 

Next, a single CMARC run is performed with a= amm and with two 
degrees of yaw angle. The results from both models are summed. The lateral-directional 
derivatives, Cyp, Qp and C n p, are then obtained directly with equations 4.22 through 4.24: 



C Y = ^ per radian 

Y * AJ3° 
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180 

71 
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7t 



per radian 
per radian 



4.22 



4.23 



4.24 



It should be noted that the CMARC wind axis is modeled with x-aft and z- 
up, vice x-forward and z-down for the flight dynamics stability axis. Care must be taken 
to reverse the signs of the appropriate coefficients to convert from the CMARC wind axis 
to the flight dynamics stability axis system. 
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Static lateral-directional stability derivatives obtained from CMARC are 
presented in Table 4.9. For comparison three other sets of data are also presented. The 
first comes from classical analysis using methods from Roskam [Ref. 12], The second 
set comes from estimates based on data recorded from flight-test sideslip maneuvers, 
published by Papageorgio in Ref. [2], The third set comes from parameter estimation, by 
Engdahl, based on dynamic flight-test maneuvers. 



METHOD 


CONFIGURATION 1 


STATIC LAT-DIR PARAMETERS 


Cyp 

(per rad) 


(per raa) 


C n [5 

(per raa) 


CMARC 
Panel Code 


Wing/Fuselage + Horz/Vert Stabs 


-0.249 


-0.063 


0.063 


Classical 

Design 2 


Wing/FuselageA/ert Tail 


-0.511 


-0.055 


0.051 


Flight Test 3 


Flying Aircraft 


-0.700 


-0.053 


0.057 


Parameter 

Estimation 4 


Flying Aircraft 


-0.987 


-0.094 


0.176 




C-172 5 


-0.310 


-0.089 


0.065 



NOTES: 1) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12], 

3) Flight test results from steady heading sideslip, from Ref. [2] 

4) Parameter estimation from flight test data by Engdahl. 

5) C-172 data from Ref. [12]. 

Table 4.9 Comparison of FROG UAV Static Lateral-Directional Stability 
Derivatives. 

CMARC produced weak results for the Cyp side force derivative. A large 
component of Cyp, approximately 60%, comes from the fuselage. Side forces are not 
modeled well by the potential flow solution. In addition, the engine pylon and pod are 
left off the model due to problems with their wakes impacting the vertical tail. They 
most likely provide a significant contribution to side force. 



54 



CMARC results for weathercock stability, C n p, and dihedral effect, Cip, 
show close correlation to the classical calculations. However, dynamic parameter 
analysis indicates a considerably higher value for both derivatives. The parameter 
estimation most likely captures the additional dynamic pressure over the vertical 
stabilizer from prop wash. As with the longitudinal solutions, an attempt should be made 
to model the propeller disk with a normal velocity distribution. 

b. Lateral-Directional Damping Derivatives 

As with the static case, the dynamic solutions are run with separate 
wing/fuselage and vertical/horizontal stabilizer models. The solutions are summed 
through superposition. Figure 4.17 shows the two angular motions selected to develop 
the lateral-directional rate-damping derivatives. The (3-dot terms are generally 
considered to be small. As a result, the sideways plunging motion test case is not run. 
Yaw rate-damping is obtained directly from the oscillating yawing motion without the 
requirement to subtract (3-dot effects. Yawing motion is conducted at a frequency of 2n 
rad/s, which equates to a reduced frequency of k=0.369 for this configuration and trim 
airspeed. 
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Roll Damping Motion 





Yaw Damping Motion 



Figure 4.17 Aircraft Motion for the Determination of Roll and Yaw damping. 
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The roll axis is de-coupled from both angle-of-attack and sideslip. 
Therefore, roll damping can be obtained directly from a constant rolling motion, as 
illustrated in Figure 4.17. For the FROG study a roll rate of 20°/sec is selected. Initially 
a lower roll rate of 5°/sec was used. However, the lower rate produced small values of 
both side force and yawing moment. A higher roll rate prevented machine resolution 
from becoming a factor in data reduction. The User Guide in Appendix A provides a 
more detailed description of the techniques used to gather the dynamic lateral-directional 
derivatives. 

The roll damping stability terms are presented in Table 4.10. Traditional 
design calculations from Roskam [Ref. 12] were completed for comparison to CMARC 
and flight-test results. In addition, C-172 data is also presented. CMARC roll damping, 
Ci p , is 50% higher than classical results, but right in line with parameter estimation. 
CMARC adverse yaw due to roll rate, C np , is approximately one fourth that obtained from 
either classical calculations or parameter estimation. Of note, CMARC does capture side 
force due to roll rate. This term is difficult to obtain from classical techniques and is 
usually assumed to be negligible. 



METHOD 


CONFIGURATION 1 


ROLL RATE TERMS 


c y P 

(per raaj 


c, P 

(per raa) 


Cnp 

(per raa) 


CMARC 
Panel Code 


Wing/Fuselage + Horz/Vert Stabs 


0.05 


-0.452 


-0.022 


Classical 

Design 2 


Wing/Fuselage/Vert Tail 


0 


-0.300 


-0.072 


Parameter 

Estimation 3 


Flying Aircraft 


0 


-0.448 


-0.108 




C-172 4 


-0.037 


-0.47 


-0.03 



NOTES: 1 ) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12]. 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref. [12]. 

Table 4.10 Comparison of FROG UAV Roll Rate Stability Derivatives. 
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The yaw rate-damping terms are presented in Table 4.11. Traditional 
design calculations from Roskam [Ref. 12] are presented for comparison to CMARC and 
flight-test results. C-172 data is also listed. CMARC does an excellent job of predicting 
yaw damping. Yaw damping, C nr , directly matches the result from parameter estimation. 
CMARC underestimates roll moment due to yaw rate, Q r , by 25 to 40 percent. And, side 
force due to yaw rate, Cy r , is overestimated by a factor of three compared to the 
parameter estimation model. No explanation is readily available to explain these errors. 



METHOD 


CONFIGURATION 1 


YAW RATE TERMS 


C Yr 

(per raa) 


c lr 

(per raa) 


Cnr 

(per raa) 


CMARC 
Panel Code 


Wing/Fuselage + Horz/Vert Stabs 


0.337 


0.121 


-0.121 


Classical 

Design 2 


Wing/Fuselage/Vert Tail 


0.140 


0.168 


-0.076 


Parameter 

Estimation 3 


Flying Aircraft 


0.110 


0.208 


-0.121 




C-172 4 


0.210 


0.096 


-0.099 



NOTES: 1) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12]. 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref [12]. 

Table 4.11 Comparison of FROG UAV Yaw Rate Stability Derivatives. 

c. Lateral-Directional Control-Power Derivatives 



Aileron and rudder control-power derivatives are presented in Tables 4.12 
and 4.13. Classical design calculations are performed using equations 4.25 through 4.30 
from Roskam [Ref. 15]: 






4.25 
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4.26 



r =r 

^'sr ^r Sr 



(Z v cos a -/ v sin a) 



C, 



”Sr 



^ (l v cos a -Z v sin a) 



4.27 



C Y& ~ small 



4.28 



=K*C l *C i 



^Sa 



4.29 



C ‘Sa ~ 



Sal 



c, 

l a 

2 71 



r c \ 

K j 



4.30 



where Cysf, oc 6 , k and (C’lga/K) are empirical values from charts in Roskam [Ref. 15]. 
Additionally, flight-test parameter analysis and C-172 data are provided for comparison. 



METHOD 


CONFIGURATION 1 


AILERO 


CONTROL POWER 


CySa 

(per raaj 


C|8a 

(per raaj 


^n5a 
(per raaj 


CMARC 
Panel Code 


Blended Wing-Fuselage/Horz/Vert Tails 


-0.021 


0.194 


-0.0121 


Classical 

Design 2 


Wing/FuselageA/ert Tail 


0 


0.213 


-0.0236 


Parameter 

Estimation 3 


Flying Aircraft 


0 


0.239 


-0.0261 




C-172 4 


0 


0.178 


-0.053 



NOTES: 1) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12], 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref. [12], 

Table 4.12 Comparison of FROG UAV Aileron Control-Power Derivatives. 
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In general, CMARC provides aileron control-power estimates in the 
correct direction and same order of magnitude as both the classical and parameter 
estimation techniques. Aileron roll control-power shows excellent correlation to both 
methods. However, CMARC under-estimates adverse yaw due to aileron deflection, 
Cn8a, by approximately 50%. Yawing moment due to roll rate, C np , is also 
underestimated by a similar margin. An investigation into the source of CMARC 
inaccuracies in the modeling of adverse yaw due to aileron deflection and roll rate is 
warranted. 

The rudder control-power derivatives obtained with CMARC show good 
correlation to classical design techniques. Side force due to rudder deflection, Cygr. is 
15% greater than classical calculations but closely matches flight-test results. For both 
C n 5 r and Cisr, CMARC estimates fall between classical calculations and parameter 
estimation results. The stronger rudder control-power observed in the parameter 
estimation model is most likely due to the capturing of prop wash effects. 



METHOD 


CONFIGURATION 1 


RUDDER CONTROL POWER 


Cy6r 

(per raa) 


C|5r 

(per raa) 


C n 5r 
(per raa) 


CMARC 
Panel Code 


Blended Wing-Fuselage/HorzA/ert Tails 


0.0928 


0.0040 


-0.0453 


Classical 

Design 2 


Wing/Fuselage/Vert Tail 


0.081 


0.0056 


-0.0341 


Parameter 

Estimation 3 


Flying Aircraft 


0.093 


0.0004 


-0.0785 




C-172 4 


0.187 


0.015 


-0.0657 



NOTES: 1) CG X =34.5% M.A.C. / CG Z =8.6” from bottom of fuselage. 

2) Classical design calculations by Roskam's methods, after Ref. [12]. 

3) Parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref [12]. 

Table 4.13 Comparison of FROG UAV Rudder Control-Power Derivatives. 



59 



3. 



Summary of CMARC Stability Derivative Analysis 



Table 4.14 lists the complete FROG aerodynamic model obtained from CMARC 
for the trim condition. In summary, CMARC produces reasonably accurate stability 
derivatives for an initial aerodynamic model. The greatest difficulty is encountered 
modeling side force due to sideslip and roll rate. The potential flow solution from 
CMARC fails to adequately capture the side force on the fuselage. Improved fidelity 
might be obtained by modeling fuselage flow separation with wake separation lines. 
CMARC also underestimates yawing moment due to roll rate and aileron deflection by 
large margins. 



LONGITUDINAL 


LATERAL-DIRECTIONAL 


Derivative 


Value 


Derivative 


Value 


c L 


0.4295 


C YP ‘ 


-0.2493 


C D 


0.065 


C|p 


-0.0630 


c La 


4.845 


Cnp 


0.6300 


Cd<x 


0.2664 


Cyp 


0.0488 


^Ma 


-0.4126 


c, P 


-0.4514 


^Ladot 


1.420 


Cnp 


-0.0220 


^Madot 


-6.264 


C Yr 


0.3370 


c Lq 


6.862 


Or 


0.1210 


^Dq 


0 


c„ r 


-0.1210 


^Mq 


-11.78 


Cysr 


0.0928 


^LSe 


0.4378 


C|5r 


0.0040 


^D6e 


0.0092 


C n £r 


-0.0453 


CfvISe 


-1.199 


Cv8a 


-0.0206 




C|Sa 


0.1943 


Cn5a 


-0.0121 



Table 4.14 Summary of CMARC Stability Derivatives for the NPS FROG UAV. 
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F. 



COMPARISON OF DYNAMIC AERODYNAMIC MODELS 



This section will discuss the dynamic response of the FROG UAV using the 
CMARC generated aerodynamic model. Specifically, modal frequency and damping are 
obtained through eigenvalue analysis. Elevator, aileron and rudder control response is 
found using a linearized state-equation model. Results are compared to the classical 
design and parameter estimation models. 

Dimensional stability derivatives are placed into linearized 4x4 models as 
outlined by Schmidt in Ref. [14]. The short-period and long-period modes are obtained 
using the MATLAB eigenvalue decomposition routine. In a similar manner, the lateral- 
directional plant matrix yields the roll, spiral and Dutch-roll modes. Dynamic response 
due to step and doublet control inputs is obtained by assembling the complete linear 
system. This paper uses control deflection sign convention consistent with Figure 4.9 
from Ref. [14]. The MATLAB "lsim" command produces the time based output of the 
linear system. MATLAB is used to plot the dynamic response. 



1. Longitudinal Dynamics 



a. Longitudinal Dynamic Modes 

The longitudinal response of an aircraft can be reduced to a series of four, 
first-order differential equations. They are typically written in matrix form based on the 
state variables u/V, a, q and 0. The linearized state-equations, in matrix form as 
developed by Schmidt in Ref. [14], are listed in equations 4.31 to 4.37: 

{*M4*}+{£)£, 4.31 

{x}=[u!V a q 0] T 4.32 



M=L!'k] 
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4.34 



vx u 


X a 


0 


-gcos© 0 


vz u 


2. 


(v + z,) 


- g sin 0 O 


VM U 


M a 


M, 


0 


0 


0 


1 


0 



4.35 



&}=[** Z„ Of 



4.36 




V 

0 

0 

0 



0 0 0 

(V-Z a ) 0 0 

-M & 1 0 

0 0 1 



Inertial matrix 



4.37 



The non-dimensional stability derivatives generated by CMARC are 
converted to dimensional derivatives through the transformations listed in Table 4.15. 
Dimensional derivatives are more convenient because they will lead to time histories 
being expressed in seconds and frequency in rad/s. The MATLAB script "froguav.m" in 
Appendix D is used to automate these calculations. 

The linearized, longitudinal, state-equation is given by equation 4.31. 
Note that the plant and control matrices are obtained by pre-multiplying by [I n ] '• The 
eigenvalues of the plant matrix, [A], yield the longitudinal dynamic modes. The 
eigenvalues are obtained from the MATLAB script "froguav.m" listed in Appendix D. 
The results are presented in Table 4.16. 
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Term 


Description 


Units 


X u 


-■5|(2C d + M^5) 

mV 3M 


s -1 


X Q 


QS, r 3 C d 
~ (Cl a* > 


ft-s -2 


x . 


QS / c \ «9 C d 


ft-s -1 


a q 


m \2 V/ <9(ac/2 V) 


x 


Q S / c \ 5Cd 


ft-s -1 




m V2vi 5(qc/2 V) 


x 6 


QS 5Cd 


ft-s -2 




m 5<5 


Zu 


Q s / 9r 

mV (2Cl+M 3M» 


s -1 


Zo 


-^(C D+ ^) 


ft-s -2 


7 • 


_ QS f c \ 5C L 


ft-s -1 


CJQ 


m V2V/ (9(orc/2 V) 


7 


QS / c \ 5C L 


ft-s -1 




m V2V/ <9(qc/2 V) 


Z* 


qs ac L 


ft-s -2 




m 5(5 


M u 


Q S c 5C m 
I y V 5M 


ft-s -1 


M 


QSc 5C m 


o 




I y da 


s 


M a 


QSC / c \ 5C m 


Q -1 


Iy V2 V/ 5(dc/2 V) 


s 


M q 


QSC / c \ 5C m 


C -1 


Iy \2V/ 5(qc/2V) 


s 


m 4 


QSc dC m 

Iv dS 


__ o 
S ~ 



Table 4.15 Dimensional Longitudinal Stability Derivatives from Ref. [14]. 
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METHOD 


CONFIGURATION 1 


LONGITUDINAL DYNAMIC MODES 


Short- 


Period 


Phugoid 


C0 n 

(raa/s) 


C 


co n 

(raa/s) 


c 


CMARC 
Panel Code 


Wing/Fuselage/Horiz Tail 


5.15 


0.909 


0.358 


0.136 


Classical 

Design 2 


Wing/Horiz Tail - 5s/8a=0.40 


5.90 


0.734 


0.407 


0.110 


Parameter 

Estimation 3 


Flying Aircraft 


4.67 


0.770 


0.397 


0.101 




C-172 


6.027 


0.685 


0.181 


0.116 



NOTES: 1) CG X =34.5% M.A.C. / CG 2 =8.6” from bottom of fuselage. 

2) Classical design calculations, after Ref. [12]. 

3) Classical model from Papageorgio [Ref. 2] as modified through 
parameter estimation from flight test data by Engdahl. 

4) C-172 data from Ref. [12]. 

Table 4.16 Comparison of FROG UAV Dynamic Longitudinal Modes. 

Table 4.16 summarizes the longitudinal dynamic modes. The CMARC 
aerodynamic model provides a better estimate of short-period frequency than the classical 
design technique. The CMARC short-period frequency is 10% faster than that observed 
from parameter analysis and damping is 21% higher. The high damping in all three 
models will produce nearly deadbeat results. The CMARC phugoid, or long period, is 
10% slower than observed from flight-test with 40% more damping. Although, it must 
be pointed out that it is difficult to accurately capture the phugoid mode in flight-test. 
This is due to difficulty the external pilot has in accurately maintaining wings level trim 
during the extended period required for capturing phugoid data. Overall, the CMARC 
aerodynamic model produces satisfactory short-period and phugoid modal data for the 
development of closed-loop flight controls. Of note, CMARC provides a more accurate 
prediction of the important short-period natural frequency. The larger error in phugoid 
frequency is less important to the design of a closed loop controller. 
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b. Longitudinal Dynamic Response to Control Input 

Following dynamic mode analysis, the response to an elevator step-input 
and a doublet-input is modeled using MATLAB. Dynamic response for the CMARC, 
classical and parameter analysis models are overlaid for comparison. The linear system 
is set up using the MATLAB script "ffoguav.m" in Appendix D. The "lsim" command 
outputs the time-based response from a time-based control input vector. A -2° (TEU) 
elevator step-input is selected to keep the response in the linear region. A 5° elevator 
doublet with a 1 .2 second period is used to excite the short-period mode. The selected 
doublet period closely matches the short-period mode. Dynamic response is displayed in 
Figures 4. 1 8 and 4.19. 

The FROG UAV response to a -2° (TEU) step elevator input is displayed 
in Figure 4.1 8. All three aerodynamic models are presented for comparison. The model 
adjusted by parameter analysis shows a much larger response than either the CMARC or 
classical design models. The CMARC model produces a final pitch angle change after 4 
seconds of 35° versus 45° for the flight-tested model. On closer review, it is noted the 
CmSe —1.199 from CMARC is considerably less than the C m s e = -1.621 from flight-test. 
It should be noted that the parameter estimation model compensates for real world factors 
including; air load distortions, sensor measurement errors and prop wash. In other words, 
the CMARC model assumes uniform displacement and perfect sensors while the flight- 
test model is empirically fit to a measured response. 

The response to a 5° elevator doublet is displayed in Figure 4.19. All 
three models show a similar frequency and a high degree of damping. As expected, the 
parameter estimation model produces a larger response. The magnitude of the CMARC 
response is approximately two thirds of that observed from the flight-test model. Again, 
this is due to the much larger value of C m s e obtained from empirically fitting the 
parameter estimation model to observed aircraft response. Modeling prop wash with 
CMARC should produce improved results. 

In summary, the CMARC aerodynamic model showed satisfactory 
longitudinal control response in frequency and damping. It is recommended that C m 5 e be 
adjusted to increase elevator response. 
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4 



LONGITUDINAL RESPONSE TO A 2 DEGREE ELEVATOR STEP INPUT 
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Figure 4.18 FROG UAV Dynamic Response to a -2° (TEU) Elevator Step Input 



66 



Theta (deg) q (deg/s) A.O.A (deg) Elevator (deg) 



LONGITUDINAL RESPONSE TO A 5 DEGREE ELEVATOR DOUBLET 
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Figure 4.19 FROG UAV Dynamic Response to a 5° Elevator Doublet 
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2 . 



Lateral-Directional Dynamics 



a. Lateral-Directional Dynamic Modes 

The lateral-directional response of an aircraft can also be reduced to a 
series of four first-order differential equations. They are typically written in matrix form 
based on the state-variables p, p, 4*, and r. The linearized state-equations in matrix form 
as developed by Schmidt in Ref. [14] are listed in equations 4.38 to 4.44: 






W=L b p 0 r Y 

M =[/„]-' kl 
tsH/.rw 




^p 

0 

N„ 



Y P 

L P 

1 

N , 



g cos© 0 (Y r + v) 
0 L r 

0 0 

0 N r 





0 N* 
0 N Sr 




V 

0 

0 

0 



0 0 0 
(F-Zj 0 0 
-M, 1 0 

0 0 1 



Inertial matrix 



4.38 

4.39 

4.40 

4.41 

4.42 

4.43 



4.44 
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As with the longitudinal axis, the stability derivatives generated by 
CMARC are dimensionless. However, it is more convenient to use dimensional 
derivatives. They will lead to time histories being expressed in seconds and frequency in 
rad/s. The lateral-directional stability derivatives are turned into dimensional derivatives 
through the transformations listed in Table 4.17 from Ref. [14]. The MATLAB script 
"froguav.m" in Appendix D is used to automate these calculations. 




Table 4.17 Dimensional Lateral-Directional Stability Derivatives from Ref. [14] 
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The linearized, lateral-directional state-equation is given by equation 4.38. 
Note that the plant and control matrices are obtained by pre-multiplying by [I n ]’'- The 
dimensional derivatives from Table 4.16 populate the linearized 4x4 plant matrix and 
control matrices. The plant matrix, [A], is used to obtain the lateral-directional dynamic 
modes. The eigenvalues are obtained from MATLAB using the "ffoguav.m" script in 
Appendix D. The results are presented in Table 4. 1 8. 



METHOD 


CONFIGURATION 1 


LATERAL-DIRECTIONAL DYNAMIC MODES 


Dutc 


h-Roll 


Roll 


Spiral 


(0 n 

(raa/s) 


c 


(0 n 

(raa/sj 


(0 n 

(raa/sj 


CMARC 
Panel Code 


Wing/Fuselage Group 
plus HorizA/ert Stab Group 


2.54 


0.13 


-3.80 


0.000 


Classical 

Design 2 


Wing/Fuselage/Horiz Tail 


2.48 


0.09 


-2.83 


0.065 


Parameter 

Estimation 3 


Flying Aircraft 


4.27 


0.14 


-3.97 


0.090 




C-172 4 


3.38 


0.20 


12.43 


0.01 



NOTES: 1) CG X =34.5% M.A.C. / CG Z =8.6" from bottom of fuselage. 

2) Classical design calculations after Ref. [12], 

3) Classical model from Papageorgio [Ref. 2] as modified through 
parameter estimation from flight test data by Engdahl. 

4) Cl 72 data from Ref. [12], 

Table 4.18 Comparison of FROG UAV Lateral-Directional Dynamic Modes. 

The lateral-directional dynamic modes are summarized in Table 4.18. The 
CMARC aerodynamic model provides significantly better estimates of frequency and 
damping than the classical design technique for both the Dutch-roll and roll modes. The 
CMARC Dutch-roll natural frequency, although better than the classical response, is 40% 
lower than predicted by parameter estimation. Adjustments should be made to the 
CMARC aerodynamic model to improve Dutch-roll response. 
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The primary contributor to the Dutch-roll frequency is weathercock 
stability, C n p. A review of Table 4.9 shows that C n p obtained from CMARC is 
considerably less than the value obtained through parameter estimation. One potential 
source of error is the lack of modeling the increased dynamic pressure due to prop wash. 
The FROG model should be re-worked to include a prop disk to investigate the ability of 
CMARC to capture prop wash effects. 

The CMARC model predicts a neutral spiral mode. This is most likely 
due to a combination of equal and opposite Q r and C nr ratios. When the Ci r /C nr ratio is 
changed to -2 (approximately that of the analytical and parameter estimation models), the 
spiral mode comes out close to the analytical model. The weak roll due to yaw-rate, Q r , 
in the CMARC solution seems to be the main source of error in the spiral mode. The 
value will need to be modified to produce an acceptable spiral response. 

b. Lateral-Directional Dynamic Response to Control Input 

Following dynamic-mode analysis, the response to aileron and rudder 
inputs is modeled using MATLAB. Dynamic response for the CMARC, classical and 
parameter-analysis models are overlaid for comparison. The linear system is set up using 
the MATLAB script "ffoguav.m" in Appendix D. The "lsim" command outputs the time- 
based response from a time-based input control vector. 

The lateral response to a +2° (right wing down) step aileron deflection is 
shown in Figure 4.20. CMARC steady-state roll-rate is less than that predicted by the 
other two models. This is expected considering the CMARC model has the lowest ratio 
of aileron control-power to roll-damping. In addition, the side-force terms, Cy p and Cysa, 
modeled by CMARC both fight roll-rate. They are not included in the other models. 
Perhaps they should be assumed small and set to zero. CMARC does a significantly 
better job of modeling the Dutch-roll response. The roll-rate and sideslip traces from the 
classical design model clearly show excess excitation of the Dutch -roll mode. The 
CMARC model shows a better correlation to the parameter estimation model for Dutch- 
roll amplitude and damping. As predicted by the eigenvalues, Dutch-roll frequency is 
slower than the parameter estimation model. 

The lateral -directional response to a 5° aileron doublet is shown in Figure 
4.21. The 1.5 second period clearly excites the Dutch-roll mode in all three models. As 
with the step-input response, the CMARC model provides a significantly better match to 
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Dutch-roll excitation than the classical design model. Amplitude and damping show a 
close correlation to the parameter estimation model. A lower Dutch-roll frequency is 
evident in the CMARC sideslip trace. 

The lateral-directional response to a +2° (nose left) rudder deflection is 
shown in Figure 4.22. The CMARC model demonstrates a close match to the classical 
design model for Dutch-roll frequency and damping. CMARC shows higher sideslip 
excitation than the classical design model. A slower Dutch-roll frequency is evident in 
both the CMARC and classical design models. 

The lateral-directional response to a 5° rudder doublet deflection is shown 
in Figure 4.23. Again, the 1.5 second doublet excites the Dutch-roll mode in all three 
models. CMARC provides a similar response to the classical design model. As 
expected, Dutch-roll natural frequency is about 40% slower than the parameter estimation 
model with similar damping. 

In summary, the CMARC lateral-directional model provides FROG 
dynamics that are similar to the classical design calculations. However, current CMARC 
lateral-directional dynamics are not adequate for closed-loop controller or autopilot 
design. Minor adjustments will be required to better match observed flight 
characteristics. The lateral-directional model should be modified to include a propeller 
disk. Higher dynamic pressure over the tail surfaces should provide a stronger C n p 
derivative and result in faster Dutch-roll frequency. In addition, the FROG model should 
be adjusted to provide a higher aileron control-power to roll-damping ratio. This 
modification will improve roll-rate response. 
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LATERAL-DIRECTIONAL RESPONSE TO A 2 DEGREE AILERON STEP INPUT 
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Figure 4.20 FROG UAV Dynamic Response to a +2° (RWD) Aileron Step Input 
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LATERAL-DIRECTIONAL RESPONSE TO A S DEGREE AILERON DOUBLET 







Figure 4.21 FROG UAV Dynamic Response to a 5° Aileron Doublet. 



74 



4 



LATERAL-DIRECTIONAL RESPONSE TO A 2 DEGREE RUDDER STEP INPUT 
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Figure 4.22 FROG UAV Dynamic Response to a +2° (TEL) Rudder Step Input 
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LATERAL-DIRECTIONAL RESPONSE TO A 5 DEGREE RUDDER DOUBLET 
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Figure 4.23 FROG UAV Dynamic Response to a 5° Rudder Doublet. 



76 



V. CONCLUSIONS AND RECOMMENDATIONS 



CMAR.C is a DOS personal-computer hosted panel-code adopted from the NASA 
Ames PMARC code. AeroLogic, Inc., created CMARC by converting PMAR.C 
FORTRAN 77 source code into the C language. Significant memory management and 
command-line enhancements were also added. CMARC solves for inviscid, 
incompressible flow over complex three-dimensional bodies. Emphasis in this study is 
placed on expanding the use of CMARC to develop a full aerodynamic model of the 
Naval Postgraduate School FROG UAV. The CMARC stability derivatives are 
compared to models derived from classical design calculations and parameter estimation. 
In addition, pitot-static and angle-of-attack sensor position corrections are obtained 
through CMARC analysis. 

The LOFTSMAN and POSTMARC portions of the Personal Simulation Works 
software suite are used exclusively for the pre-process modeling and post-process 
visualization of CMARC files. The LOFTSMAN capability to automatically format and 
generate CMARC input patches is an enhancing characteristic. Functionality should be 
added to allow the modeling of wing tip ribs that are not parallel to the aircraft centerline. 

POSTMARC is an excellent tool for visualizing CMARC output files. The 
capability to create streamlines and perform boundary- layer calculations external to 
CMARC is extremely useful. However, much time could be saved if POSTMARC 
maintained previous settings and selections following translations, rotations and re- 
scaling. Additionally, a capability to overlay multiple data types is desired. 

CMARC off-body flow-field analysis is useful for both static source and angle-of- 
attack sensor position corrections. In-flight measurements may be corrected using look- 
up tables or through curve-fits of CMARC-derived data. Flight testing is recommended 
for validation of sensor corrections obtained from the CMARC off-body analysis. 

For the static longitudinal analysis, CMARC produces accurate values for a mm 
and C Ma and a slightly high value of C La . Elevator control-power (C M5e ) from CMARC is 
considerably weaker than the value obtained from parameter estimation. CMARC- 
derived pitch-damping is stronger than the pitch-damping from both classical design 
calculations and parameter estimation. 

Longitudinal dynamic-mode analysis shows an acceptable match for short-period 
and long-period frequency and damping. As expected, a lower dynamic response to 
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elevator control inputs is observed when compared to the parameter estimation model. 
However, in all cases, the CMARC aerodynamic model demonstrates better dynamic 
response than the classical design model. A propeller disk should be added to the 
CMARC model in an attempt to capture prop-wash effects on the horizontal stabilizer 
and elevator. Increased elevator control-power would provide a better match with 
observed flight characteristics. 

In general, CMARC provides satisfactory lateral-directional derivatives. Roll- 
damping closely matches the value obtained through parameter estimation. The greatest 
difficulty encountered is in modeling side-force due to both sideslip and roll-rate. The 
potential flow solution from CMARC fails to adequately capture the side force on a 
slender-body fuselage. Additional wakes, placed along separation lines, may improve 
fuselage side-force prediction. 

Dynamic response to rudder control input shows close correlation in amplitude 
and damping to the parameter estimation model. However, both the CMARC and 
classical models of the FROG UAV produced a slower Dutch roll frequency. Steady 
state roll-rate obtained from the CMARC model is somewhat slower than either the 
classical design or parameter estimation models. 

Overall, the CMARC panel code is found to be suitable for aerodynamic 
modeling of the Naval Postgraduate School FROG UAV. CMARC-derived stability 
derivatives are sufficiently accurate for incorporation into an initial aerodynamic model. 
The CMARC aerodynamic model demonstrated better longitudinal dynamic response 
than the classical design model. Lateral-directional response is similar to that obtained 
from classical design techniques. Adjustment through comparison with flight-test data is 
still required to optimize the CMARC model. Future studies should concentrate on 
improving CMARC modeling of fuselage side force through the addition of separation 
wake lines. Additionally, the propeller disk should be modeled in an attempt to capture 
the effects of increased dynamic pressure over the horizontal and vertical tail surfaces. 
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A. 



INTRODUCTION 



1. CMARC Panel Code and Limitations 

CMARC is a PC hosted panel code. The code is useful for gathering aircraft 
stability data subject to the limitation of inviscid, potential flow solutions. Panel codes 
assume fully attached flow. Therefore, care must be taken to gather stability data in the 
linear region of the lift curve slope where separation effects are minimized. The results 
obtained from CMARC should always be spot checked with classical design calculations. 

CMARC and PMARC have been shown to produce functionally equivalent 
results. The guide describes the use of DOS batch files for automating CMARC data 
runs on a PC. PMARC, hosted on an SGI machine, will require the use of UNIX batch 
commands. 

2. Overview of Obtaining Stability Derivatives with CMARC 

An overview of the methods to obtain stability derivatives follows. The detailed 
description and data reduction methods can be found in the appropriate section of the user 
guide. 

For the static stability derivatives, the CMARC model is run at two different 
angles-of-attack and one sideslip angle. The static solutions are obtained without control 
surface deflections. The slopes of the force and moment coefficients are then taken to 
produce the C La and C ma longitudinal derivatives and the C Yp , C ip and C„ p lateral- 
directional derivatives. 

For the control power derivatives, the model is run at the trim condition with 
successive elevator, aileron and rudder control deflections. The difference between the 
results with and without control surface deflections yield the C LSe , C M5e and C DSe 
longitudinal and the C Y5r , C )5r , C n8r , C 15a and C n5a lateral-directional control power 
derivatives. 

Development of the damping derivatives is not as straight forward. For the 
dynamic derivatives, motion is enabled around the center of gravity. With motion turned 
on, it is extremely important to ensure that the global origin and C.G. are co-located. If 
they are not, a rotation around a specific axis will also create a translation at the C.G. 
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For the longitudinal damping derivatives, the model is run with oscillating vertical 
plunging motion to obtain the C L and C M a-dot terms. The lift and pitching moment 
coefficients are broken into real (in phase with AOA) and imaginary (out of phase with 
AOA) components. The imaginary components are due to a-dot effects. Next, the model 
is run with oscillating pitch motion to obtain the combined a-dot and pitch rate terms. 
Subtracting the a-dot influence obtained from the plunging motion isolates the pitch rate 
damping term from the pitch motion. 

For the lateral-directional analysis, the P-dot terms are generally negligible. This 
allows the model to be run with just oscillating roll and yaw motion. As with the 
longitudinal test case, the imaginary or out of phase component yields the combined P- 
dot and rate damping data. With the p-dot terms assumed negligible, the oscillating 
motion yields the C Y , C„ and C n roll and yaw rate terms directly. 

3. Coefficient Normalization and Stability Axis System 

CMARC contains built-in functionality to integrate forces and moments over the 
surface of a body. Forces and moments are automatically normalized into non- 
dimensional coefficients based on the mean aerodynamic chord, reference wing area, 
semi-span and center of gravity location in the CMARC BINP9 input line. CMARC 
outputs coefficients in both the wind and body axes. Of note, CMARC uses the semi- 
span to normalize rolling and yawing moment coefficients. Most texts on stability and 
control, including Roskam's "Aircraft Flight Dynamics" and Etkin's "Dynamics of 
Flight," normalize rolling and yawing moments by span. This user guide will normalize 
roll and yaw moments by span. Table A.l summarizes the factors for normalizing 
moments and angular rates. All rolling and yawing moment coefficients presented in this 
guide have been normalized with span by dividing the CMARC output by a factor of two. 
Table A.l also indicates the characteristic time, t*, employed for angle rate data 
reduction. 

In addition to differences in normalizing moments, CMARC uses the typical CFD 
axes system shown in Figure A.l. This user guide will perform all calculations in the 
stability axes system as illustrated in Figure A.l. The sign of CMARC roll and yaw 
moments need to be reversed. The direction for positive control deflections is also shown 
in Figure A.l. All control surfaces are patched with positive deflections using the 
convention in Figure A.l. 
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Note: 1 ) CMARC normalizes roll and yaw coefficients with b/2. 



Table A.l Normalized Moment and Rate Equations. 

B. MODEL VALIDATION 

1. Initial Solution at the Trim Condition 

Once a CMARC model is built up, a reference solution is obtained at the desired 
trim condition. Initial hand calculations or flight test data can be used to select the 
estimated cruise angle-of-attack. Use POSTMARC to visualize the reference solution. 
Check to make sure the pressure distribution is consistent with expected results. 
Stagnation zones (Cp=l) should be evident on the leading and trailing edges of the flight 
surfaces and aerodynamic bodies. Lower pressure should be observed on the top of 
lifting surfaces and tip loss should be observed near the wing tips. Wakes should trail 
from the expected surfaces. 

If abnormal pressure distribution is observed, check for inverted panel orientation. 
This is accomplished by changing the color for positive and negative panel orientations 
using the "View » Color" pull-down menu. The default color is white for both 
orientations. Change the inside to a contrasting color such as red. All panels should face 
in the positive out direction. If a patch has panels oriented inward, change the orientation 
of the patch with the "IREV" term in the initial patch definition. 

Check that the model is symmetric around the lateral and directional axes with 
sideslip set to zero. Residual rolling or yawing moment without a control deflection is an 
indication of an asymmetric model. 
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,2 




CMARC Axis System 




Stability Axis System 



Figure A.l CMARC CFD Axis System Compared To Stability Axis System. From 
References (3] and [14] Respectively. 



86 



Pitching moment should be checked for a reasonable value. Don't forget to verify 
that the origin of the global coordinate system is set to the center-of-gravity (C.G.). It is 
important to remember that the damping derivatives will be obtained by inducing motion 
around the center-of-gravity. If the C.G. and origin are not collocated, angular motion 
will generate undesired translation at the C.G. 

Next, the lift coefficient from the initial solution is checked to make sure it is in 
the "ballpark." If a major discrepancy is noted, the input file will have to be checked 
closely for error. Be sure to check the normalizing factors such as wing area, reference 
chord and semi-span. If every thing looks good, try increasing or decreasing the angle- 
of-attack to see if a slight adjustment will bring the lift coefficient in line with the hand 
calculation. Once the model exhibits symmetry with a reasonable lift distribution, you're 
ready to start gathering data. 
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c. 



LONGITUDINAL STABILITY DERIVATIVES 



This section will describe CMARC methods for the development of longitudinal 
stability derivatives to include the static, rate damping and control power derivatives. As 
stated earlier, the results obtained from CMARC should be checked against classical 
design calculations. The potential flow analysis performed by CMARC does not provide 
accurate viscous drag values. Total drag can be estimated from the flight test power-off 
glide ratio or cruise thrust required. If this data is unavailable, then empirical data from 
publications such as Hoemer's "Fluid Dynamic Drag" can be used to estimate total drag. 

For the longitudinal analysis, only half the model is analyzed. This cuts the 
model size in half, resulting in much quicker solutions. The symmetric calculation mode 
is selected by setting both RSYM=0.0 and EPATSYM=0 in the CMARC input file. 
Remove the vertical tail patch if it interferes with wing or fuselage wakes. Figure A.2 
shows the FROG model configuration used to find the longitudinal derivatives. A rigid 
wake is selected that is continuous from wingtip to wingtip. As will be seen in a later 
section, the vertical tail is added for the lateral-directional solutions. This requires the use 
of a modified wing wake to avoid the vertical tail. 





Blended Wing and Fuselage with Horizontal Tail 



Figure A.2 FROG UAV Model for Obtaining Longitudinal Stability Derivatives. 
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1. 



Static Longitudinal Stability Derivatives 



a . Static Longitudinal Stability Equations 



Three basic longitudinal stability derivatives can be measured with just 
two runs of the CMARC model. The model is first analyzed at an angle-of-attack 
corresponding to the estimated trim condition. For the FROG case, a =0° is selected for 
the first run based on hand calculations. A second CMARC run is conducted with angle- 
of-attack incremented by one or two degrees. C L , C D and C m are then extracted manually 
from the data files. The slope of C L , C D and C m versus angle-of-attack provides the C La , 
C Da and C ma longitudinal derivatives. Additional solutions should be obtained to check 
for lift slope linearity. In addition, a tIim is calculated from the lift curve slope and trim lift 
coefficient. Equations A. 1 through A.4 are used for these calculations. 



Cl — 

L a 



(Cl 2 -oj 180 



C D a = 



r = 

^m r 



(«2 - 


-«l) 


71 


[C0 2 


~ C D\) 


,180 


(cc 2 


-«i) 


71 


{Cm 2 


-c*. J 


1,180 


(« 2 


-a x ) 


71 



per radian 



A.l 



A.2 



A.3 



cr 



trim 



= a° i + 



1 ' C krim - C £l)*180 






* degrees 

71 



A.4 



b. Sample Static Longitudinal Stability Data Reduction 

Table A.2 presents CMARC solutions for the FROG UAV at two angles- 
of-attack. Sample calculations for obtaining the static derivatives are demonstrated 
below. The calculations are easily implemented in a spreadsheet or with a MATLAB 
script. 
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RUN# 


AOA 


CL 


CD 


CY 


C m 


C n 


C 1 


1 


0° 


0.4300 


0.0277 


0.0000 


0.0444 


0.0000 


0.0000 


2 


2° 


0.5991 


0.0370 


0.0000 


0.0300 


0.0000 


0.0000 



Table A.2 FROG UAV Static Longitudinal Stability Data. 



Sample Calculations: 
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c. Sample Total Drag Coefficient Calculation 

The final static longitudinal parameter required is total aircraft drag. C D 
plays an important role in long period aircraft dynamics. Unfortunately, potential flow 
panel codes such as CMARC do not provide accurate total drag estimates. They can 
provide good induced drag predictions. In addition, if equipped with a boundary layer 
code like that contained in CMARC, they can provide integrated skin friction results. 
However, a large total drag contribution in the form of separation drag is not accounted 
for. Total drag estimates are made below using the two simple techniques shown in 
Equations A.9 to A. 12. The first method is based on the flight-tested glide ratio. The 
second is based on cruise power required and estimated prop efficiency. Note that the 
selected prop efficiency is relatively low due to the FROG UAV small propeller diameter, 
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